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ABSTRACT 


A  high  temporal  resolution  data  set  from  a  mooring  in  Monterey  Bay,  California  was  analyzed 
and  used  to  calculate  heat  and  momentum  fluxes  for  the  purpose  of  forcing  two  ocean  mixed  layer  models. 
The  time  frame  for  the  study  was  September.  1992.  a  period  representative  of  the  sea  breeze  circulation 
frequently  affecting  this  and  other  coastal  regions. 

The  models  used  were  that  of  Price.  Weller  &  Pinkel  ( 1986),  a  Richardson  number  based 
mixing  model,  and  Garwood  (1977').  a  model  based  on  the  turbulent  kinetic  energy  budget  within  the 
mixed  layer.  Both  models  were  analyzed  with  respect  to  their  ability  to  reproduce  the  observed  diurnal 
variation  of  the  temperature  and  depth  of  the  mixed  layer.  Although  the  model  predictions  agree 
reasonably  well  with  observations  in  regards  to  the  phase  of  the  diurnal  temperature  cycle,  they  were  seen 
to  underpredict  its  magnitude,  particularly  the  nocturnal  cooling.  This  lack  of  cooling  in  the  models 
relative  to  the  ocean  could  be  due  to  penetrative  convection,  non-steady  state  turbulence,  and/or  diurnal 
advection  present  in  the  ocean  but  not  in  one  or  both  models.  Additionally,  the  models  exhibited  an 
upward  temperature  trend  relative  to  the  data  which  caused  progressively  increasing  stratification.  This 
trend  was  used  to  approximate  the  magnitude  of  vertical  advective  effects. 
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I.  INTRODUCTION 


The  uppermost  layer  of  the  ocean  is  commonly  referred  to 
as  the  surface  mixed  layer,  due  to  its  relatively  uniform 
velocity  structure  and  isothermal,  isohaline  characteristics. 
The  depth  of  the  mixed  layer  is  primarily  determined  by  the 
balance  between  turbulence  generated  at  the  surface  by  the 
wind,  and  buoyant  fluxes  caused  by  surface  heating  and 
cooling.  The  wind  acts  to  generate  turbulence,  eventually 
resulting  in  entrainment  of  water  from  below  the  mixed  layer 
and  a  corresponding  cooling  and  deepening  of  the  layer. 
Conversely,  when  the  buoyant  flux  is  positive  downward,  as  is 
the  case  on  a  sunny  summer  afternoon  in  the  mid- latitudes,  the 
mixed  layer  will  shallow  and  warm.  This  occurs  due  to 
increased  stratification  of  the  uppermost  portion  of  the  water 
column,  unless  the  wind  mixing  is  strong  enough  to  overcome 
this  effect.  Thus,  in  the  one -dimensional  case,  the  balance 
of  these  forces,  shear  instability  due  to  wind,  and  buoyant 
forces  due  to  surface  heat  flux,  is  the  principal  factor  in 
determining  the  depth  and  temperature  of  the  ocean  mixed  layer 
in  a  given  region.  Other  factors  affecting  mixed  layer  depth 
and  temperature  are  the  input  of  fresh  water  from  river 
runoff,  precipitation,  evaporation,  and  shear  stress  at  the 
base  of  the  mixed  layer  due  to  the  presence  of  internal  waves. 
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Of  course,  advection  of  heat,  salt,  and  momentum  add  further 
complexity  to  this  one -dimensional  view. 

Studies  of  mixed  layer  dynamics  have  generally  involved 
long  time  scale  forcing  in  open  ocean  regions.  That  is,  they 
have  been  concerned  with  synoptic  scale  or  larger  weather 
patterns  and  monthly  or  seasonal  insolation  patterns.  This 
study  uses  much  of  the  knowledge  of  mixed  layer  dynamics  and 
air/sea  fluxes  gained  in  these  open  water  efforts  and  extends 
it  to  a  coastal  region.  Specifically,  the  area  of  interest  in 
this  paper  is  Monterey  Bay,  off  central  California.  The 
observed  diurnal  cycle  in  the  radiation  here  is  similar  to 
previous  studies,  but  the  wind  stresses  also  undergo  a 
strongly  diurnal  variation,  peaking  each  afternoon  in  what  is 
commonly  referred  to  as  a  land/sea  breeze  circulation.  During 
much  of  the  primary  period  of  study,  in  fact,  the  area  weather 
is  under  the  influence  of  a  quasi -stationary  high  pressure 
system,  which  greatly  reduces  wind  variations  on  a  synoptic 
scale.  This  leaves  the  local  land/sea  breeze  influences  and 
daily  heating  as  the  primary  controllers  of  mixed  layer 
variability,  with  the  advection  of  cooler  water  into  the  study 
area  playing  a  role  as  well  in  this  upwelling  favorable 
eastern  boundary  current  regime.  Rainfall  is  rare  during  the 
summer  months  in  the  study  area  (none  fell  during  the  period 
of  interest)  and  river  runoff  is  very  slight,  allowing  fresh 
water  input  to  be  neglected  in  the  modeling  process.  Also  not 
considered  in  the  models  are  the  shear  stresses  at  the  base  of 
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the  mixed  layer,  since  these  were  considered  to  be  small  in 
comparison  to  heat  flux  and  surface  shear. 

This  project  was  undertaken  in  an  effort  to  broaden 
knowledge  of  coastal  ocean  dynamics  as  part  of  the  Naval 
Meteorology  and  Oceanography  Command's  recent  shifting  of 
emphasis  from  the  deep  ocean  to  nearshore,  in  accordance  with 
the  CNO  doctrine  set  forth  in  From  the  Sea.  In  particular, 
this  study  will  give  insights  into  mixed  layer  behavior  in  a 
region  in  which  diurnal  wind  variability  is  a  dominant  feature 
of  the  overall  wind  stress  pattern.  In  addition  to  its 
contribution  to  the  science  of  coastal  ocean  dynamics  and 
mixed  layer  processes,  this  effort  will  have  applications  to 
naval  operations,  including  coastal  ASW,  mining,  diving,  and 
amphibious  operations,  particularly  in  areas  with  significant 
sea  breeze  signatures,  which  includes  the  coasts  of  most  of 
the  low  and  mid  latitudes.  The  area  of  study  was  selected 
because  of  the  availability  of  a  unique  data  set  with  high 
temporal  resolution,  combined  with  a  wealth  of  local  area 
knowledge  and  supplementary  data. 

Few  studies  have  examined  the  shallow  ocean  response  to 
diurnal  wind  forcing.  Rosenfeld  (1988)  presents  results  from 
the  Coastal  Ocean  Dynamics  Experiments  (CODE)  I  and  II,  which 
indicate  that  off  Point  Arena  in  northern  California  there  is 
a  significant  diurnal  variation  to  the  wind  induced  currents. 
The  CODE  region  is  also  characterized  by  a  predominantly 
diurnal  wind  variation  during  summer  months,  but  the  wind 
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accelerates  along  the  coast  to  the  south  during  the  day  in 
response  to  large  scale  heating  in  the  central  valley  of 
California  (Beardsley  et  al .  ,  1987)  .  In  contrast,  the  present 
study  area  undergoes  more  of  a  classical  onshore/offshore  wind 
variation  due  to  more  local  heating  in  the  Salinas  Valley  and 
the  winds  are  directionally  controlled  by  the  northwest  to 
southeast  orientation  of  that  valley  (Figure  1)  .  While  the 
coastal  mountains  in  the  CODE  area  are  basically  continuous, 
a  significant  break  occurs  in  the  Salinas/Monterey  area, 
allowing  the  diurnal  cross -shore  flow  to  develop  through  much 
of  the  year.  Other  unique  aspects  of  the  area  are  the 
presence  of  strong  upwelling  centers  to  the  north  and  south  of 
Monterey  Bay  and  the  Monterey  Submarine  Canyon,  which  cuts 
through  the  center  of  the  bay.  The  upwelling  that  occurs  in 
the  area  to  the  north  of  our  study  area  and  is  advected 
southward  (Rosenfeld  et  al.,  1993)  and  the  possible  upwelling 
at  the  study  site  add  to  the  complexity  of  this  coastal  region 
by  potentially  introducing  significant  horizontal  and  vertical 
advection  to  the  mixed  layer  problem. 

The  mixed  layer  of  Monterey  Bay  will  be  investigated 
through  a  careful  analysis  of  the  available  data  during  1 
through  11  September  1992  and  a  comparison  of  these  data  to 
results  produced  by  two  one -dimensional  mixed  layer  models, 
those  of  Garwood  (1977)  and  Price,  Weller  &  Pinkel  (1986)  . 
These  two  models  represent  two  very  different  approaches  to 
mixed  layer  physics.  The  heat  and  momentum  fluxes  that  drive 
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the  models  are  computed  from  the  measured  data.  The  models' 
ability  to  reproduce  the  observed  diurnal  cycle  of  mixed  layer 
temperature  and  depth  will  be  of  primary  interest,  with  a 
secondary  goal  of  determining  how  well  they  reproduce  the 
trend  over  an  eleven  day  period.  This  latter  effort  is  made 
in  an  attempt  to  quantify  the  advection  necessary  to  maintain 
the  cool  temperatures  of  the  bay  in  the  presence  of  large 
downward  heat  fluxes. 

Data  used  in  this  study  are  largely  from  mooring  Ml  (in 
the  center  of  the  bay  as  seen  in  Figure  1)  ,  owned  and  operated 
by  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI)  .  This 
mooring  collects  a  suite  of  meteorological  and  ocean 
temperature  data  at  ten  minute  intervals.  These  data  were 
provided  courtesy  of  Francisco  Chavez  of  MBARI.  Ocean 
velocities  measured  by  a  downward  looking  Acoustic  Data 
Current  Profiler  (ADCP)  at  fifteen  minute  intervals  was 
provided  by  Leslie  Rosenfeld  of  MBARI  and  the  Naval 
Postgraduate  School  (NPS) .  These  instruments  will  be 
described  in  greater  detail  in  the  next  section.  Some 
meteorological  data  were  also  obtained  from  the  NPS 
meteorological  station  at  Fritzche  Field  on  Fort  Ord  (Figure 
1)  . 

Some  basic  concepts  of  ocean  mixed  layer  physics  will  be 
presented  in  the  next  chapter,  followed  by  a  brief  account  of 
previous  one -dimensional  modeling  efforts  and  coastal  mixed 
layer  studies.  Chapter  III  is  a  presentation  of  the  key 
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features  of  the  two  models  used  in  this  study.  A  description 
of  the  data  sources  and  the  data  itself  is  contained  in 
Chapter  IV,  including  a  presentation  of  computed  fluxes,  while 
Chapter  V  is  a  description  of  the  model  sensitivities  along 
with  the  results  of  the  model  runs.  Chapter  VI  is  the 
discussion  portion  of  the  thesis,  with  detailed  analysis  of 
the  performance  of  the  models  compared  with  the  observations. 
Finally,  conclusions  and  recommendations  for  further  research 
and  model  improvement  will  be  presented  in  Chapter  VII. 
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Longitude  (-=West) 


Figure  1.  Monterey  Bay,  California  and  the  Salinas 
Valley.  Mooring  (Ml)  and  Profiler  (PRO)  Locations  are 
Shown.  Also  Shown  for  Reference  are  Monterey  (MBA)  and 
Santa  Cruz  (UCSC) .  Contour  Interval  is  200  Meters. 
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II.  BACKGROUND 


A.  MIXED  LAYER  BASICS 

It  is  now  widely  accepted  that  in  all  areas  of  the  world's 
oceans,  there  exists  a  layer  of  varying  depth  near  the  surface 
which  may  be  considered  as  a  nearly  isothermal,  isohaline  and 
isovelocity  slab.  This  allows  temperature,  salinity  and 
velocity  within  the  mixed  layer  to  be  approximated  as  bulk 
quantities,  greatly  simplifying  the  governing  equations.  A 
simple  mixed  layer  diagram  which  reflects  this  bulk  picture  is 
presented  in  Figure 
2 .  The  depth  at 
which  the  sharp 
gradients  of  these 
variables  begins  is 
referred  to  as  the 
mixed  layer  depth, 
termed  h.  The  thin 
layer  below  this 
depth  is  called  the 
entrainment.  zone, 
with  thickness  5. 

This  is  the  region 
in  which  cold, 


Figure  2.  Bulk  Mixed  Layer  Diagram. 
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dense  water  from  below  is  entrained  into  the  mixed  layer, 
lowering  the  mixed  layer  temperature. 

The  mixed  layer  depth,  and  thus  its  temperature,  is 
determined  by  the  amount  of  turbulence  available  for  mixing, 
as  discussed  in  the  introduction.  Shear  instability  at  both 
the  top  and  bottom  of  the  mixed  layer  contributes  to  this, 
although  the  surface  forcing  due  to  wind  stress  appears  to  be 
of  much  greater  significance.  As  the  wind  stress  acting  on 
the  ocean  surface  increases,  turbulence  generation  due  to 
shear  also  increases.  Heat  flux  at  the  surface  is  the  other 
key  factor.  When  the  total  flux,  a  combination  of  heat  gained 
from  solar  insolation,  Qg,  and  the  heat  lost  from  the  surface 
due  to  latent,  Qe<  sensible,  Qh,  and  net  infrared  radiation, 
Qb  is  positive  downward,  as  would  generally  be  the  case  during 
the  day  in  mid- latitude  summer,  mixed  layer  turbulence  is 
damped  as  the  stratification  increases.  Bulk  theory  assumes 
that  the  water  below  the  entrainment  zone  remains  completely 
non- turbulent .  At  night,  the  heat  loss  from  the  surface 
becomes  dominant  in  the  absence  of  downward  heat  flux  (in 
general)  and  buoyant  turbulence  generation  occurs,  enhancing 
any  wind  mixing  which  may  be  present.  This  occurs  through 
convective  overturning  of  the  water  as  its  surface  is  cooled. 
In  summary,  wind  shear  at  the  surface  always  acts  to  generate 
turbulence  which  is  available  for  mixed  layer  deepening,  while 
buoyant  fluxes  may  be  either  a  source  or  a  sink  for  the 
generation  of  turbulence. 
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B.  MIXED  LAYER  MODELING 


Treatment  of  the  mixed  layer  with  a  one-dimensional  model 
is  a  valid  approach,  in  that  vertical  variations  in  water 
properties  over  about  100  meters  typically  far  exceed  any 
horizontal  variations  over  1000  km  or  more  (Niiler  and  Kraus, 
1977) .  This  allows  horizontal  derivatives  to  be  neglected, 
greatly  simplifying  the  dynamics  needed  in  the  model.  An 
early  mixed  layer  model,  Ball  (1960),  which  was  actually 
developed  for  the  atmosphere,  studied  convective  effects  in 
the  absence  of  horizontal  motion  over  heated  ground.  This  was 
extended  to  ocean  mixed  layer  applications  by  the  pioneering 
work  of  Kraus  and  Turner  (1967) ,  whose  model  became  the  basis 
for  the  group  of  so  called  integrated  models,  which  are  based 
on  the  bulk  assumptions  mentioned  previously  and  which  include 
both  of  the  models  used  in  this  study. 

There  are  basically  three  other  categories  of  mixed  layer 
models  as  described  by  McCormick  and  Meadows  (1988)  .  The 
first  of  these  is  the  turbulence  closure  models,  such  as 
Mellor  and  Yamada  (1974,  1982)  and  Kundu  (1980) ,  which  solve 
for  Reynolds  stress  terms  in  the  turbulent  energy  equations 
through  higher  order  terms.  These  models  are  relatively 
complex  and  require  additional  assumptions  and  empirically 
defined  constants.  Next  are  the  deterministic  solutions,  such 
as  Deardorff  (1970) .  This  model  computes  the  Reynolds  terms 
directly  from  variables  which  must  be  known  with  very  fine 
spatial  and  temporal  resolution  and  is  extremely  time 
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consuming.  Finally,  there  are  the  eddy  diffusion  models, 
including  Kent  and  Pritchard  (1959),  Pacanowski  and  Philander 
(1981),  and  McCormick  and  Scavia  (1981).  These  are  based  on 
the  thermal  energy  equation  and  the  assumption  that  the 
Reynolds  terms  can  be  expressed  according  to  a  local 
relationship  between  mean  scalar  fields  and  eddy  fluxes,  as 
first  done  by  Munk  and  Anderson  (1948).  Martin  (1985) 
demonstrated,  using  a  comparison  of  the  Garwood  (1977) ,  Niiler 
(1975),  and  Mellor  and  Yamada  (1974,  1982)  models  that,  at 
ocean  weather  stations  Papa  and  November  in  the  Pacific,  the 
complicated  models  of  the  turbulence  closure  type  did  not 
perform  better  in  any  significant  way  than  the  simpler 
integrated  models.  In  fact,  the  Garwood  model  reproduced  the 
temperature  pattern  at  both  stations  better  than  the  others. 
The  deterministic  models  have  proven  too  unworkable  and  the 
eddy  diffusion  type  has  received  much  criticism.  Therefore, 
this  study  focused  on  two  models  of  the  integrated  Kraus  and 
Turner  type,  which  will  be  discussed  here  in  more  detail. 

The  Kraus  and  Turner  model  was  seen  in  several  studies  to 
have  a  problem  with  excessive  wintertime  erosion  of  the 
thermocline  and  it  was  widely  considered  that  an  improvement 
was  necessary  in  the  area  of  the  parameterization  of  viscous 
dissipation  (Martin,  1985) .  In  other  words,  by  neglecting  the 
viscous  dissipation,  the  model  retained  excessive  turbulence 
and  thus  overdeepened  the  mixed  layer.  The  models  of  Geisler 
and  Kraus  (1969),  Miropol'skiy  (1970),  Denman  (1973),  and 
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Niiler  (1975)  were  all  variations  of  Kraus  and  Turner  (1967) 
with  dissipation  included  as  a  fixed  fraction  of  wind  stress 
production  (Garwood,  1977) .  Resnyanskiy  (1975) ,  Kim  (1976) , 
and  Elsberry  et  al .  (1976)  extended  this  by  recognizing  the 
need  for  additional  dissipation  during  certain  types  of  storm 
forcing.  Finally,  Pollard,  Rhines  and  Thompson  (1973) 
deviated  from  other  Kraus  and  Turner  based  models  by  utilizing 
a  total  kinetic  energy  budget,  rather  than  a  turbulent  kinetic 
energy  budget.  Whereas  the  original  Kraus  and  Turner  type 
determines  entrainment  by  weighing  the  wind  generated 
turbulence  and  buoyant  forces  at  the  base  of  the  mixed  layer, 
this  second  type,  referred  to  as  the  dynamic  instability  type 
by  Cushman -Rois in  (1981)  ,  deepens  the  layer  when  the  mean  flow 
becomes  unstable  based  on  some  criterion.  The  two  models 
presented  in  this  study  represent  these  two  branches  of  Kraus 
and  Turner  based  integrated  mixed  layer  models.  Garwood's  is 
of  the  turbulence  budget  type,  while  Price  et  al .  (1986)  is 
based  on  mean  flow  dynamic  instability  for  its  mixing. 

C.  MIXED  LAYER  BEHAVIOR  IN  COASTAL  REGIONS 

As  mentioned  in  the  introduction,  most  of  the  above  models 
have  been  applied  principally  to  mid- ocean,  long  time  scale 
problems.  There  have  been  some  studies,  however,  in  which 
coastal  influences  and  diurnal  scale  mixed  layer  behavior  have 
been  investigated.  Price  et  al.  (1986)  looked  at  the  diurnal 
cycle  in  the  temperature  and  currents  using  high  temporal 
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resolution  data  from  R/P  Flip.  They  defined  a  "trapping 
depth"  as  the  mean  depth  value  of  a  temperature  anomaly 
profile,  where  the  temperature  anomaly  is  the  difference 
between  the  temperature  at  a  given  depth  and  some  reference 
temperature.  They  observed  that  the  key  to  understanding  the 
diurnal  cycle  of  mixed  layer  heat  content  was  to  determine  how 
this  depth  in  responds  to  surface  wind  stress  and  stabilizing 
surface  heat  flux.  Their  mixed  layer  model  will  be  presented 
in  the  next  chapter. 

The  Coastal  Upwelling  Ecosystems  Analysis  (CUEA)  program 
in  the  1970 's  resulted  in  several  studies  in  which  a  diurnally 
varying  wind  was  observed.  Like  the  winds  of  the  northern 
California  CODE  region,  the  winds  off  Oregon  and  Peru  were 
found  by  Halpern  (1974)  ,  Burt  et  al.  (1973,  1974) ,  and  Johnson 
(1975)  to  vary  mainly  in  the  along  shelf  direction,  whereas 
the  winds  off  Africa  were  more  similar  to  the  Monterey  Bay 
region  in  that  they  exhibited  more  of  a  cross  shore  sea  breeze 
type  circulation  (Halpern,  1977) .  All  three  regions  exhibited 
a  clockwise  rotation  of  the  diurnal  winds  (Rosenfeld,  1988) . 
Off  Africa,  a  diurnal  cross -shelf  wind  magnitude  of  3.5  m/s 
was  observed  out  to  32  km  from  shore  (Halpern,  1977) .  These 
studies  also  observed  that  the  surface  currents  rotated 
clockwise,  with  the  winds.  Halpern  (1977)  concluded  that 
these  diurnal  currents,  since  they  were  only  in  the  upper 
water  column,  might  be  generated  by  the  diurnal  period  wind, 
but  no  direct  connection  had  been  established  (Rosenfeld, 
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1988) .  Rosenfeld  (1988)  studied  the  surface  currents  in  the 
CODE  area  and  modeled  mixed  layer  behavior  using  the  model  of 
Price  et  al .  (1986).  She  observed  diurnal  currents  2  to  3 
times  their  average  value  during  times  of  strong  upwelling 
favorable  winds.  Again,  they  were  found  to  be  clockwise 
rotating  and  surf ace- intensified,  with  magnitudes  up  to  20 
cm/s.  A  strong  correlation  with  the  diurnal  winds  was 
observed. 

These  previous  studies  have  not  thoroughly  investigated 
the  diurnal  mixed  layer  thermal  structure.  The  Rosenfeld 
(1988)  work  focused  primarily  on  describing  upper  ocean 
currents  at  the  diurnal  period.  Price  et  al .  (1986)  looked  at 
the  diurnal  temperature  cycle,  but  not  in  a  region  with 
diurnally  varying  winds.  The  present  study  is  focused  on  this 
thermal  structure  in  a  sea  breeze  influenced  region.  Currents 
are  not  specifically  analyzed  in  this  study,  as  a  one- 
dimensional  model  is  not  generally  as  useful  for  this  purpose 
as  it  is  for  predicting  temperature  variations  (Niiler  and 
Kraus,  1977) .  Also,  companion  studies  to  this  one  (Foster, 
1993  and  Petruncio,  1993)  investigate  thoroughly  the  diurnal 
wind-driven  and  tidal  currents,  respectively,  within  Monterey 
Bay. 
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III.  MODEL  DESCRIPTIONS 


Both  the  Garwood  (1977)  and  the  Price  et  al .  (1986)  models 
(hereinafter  referred  to  as  Garwood  and  FWP)  are  one¬ 
dimensional  ocean  mixed  layer  models.  That  is,  they  use  only 
vertical  heat  and  momentum  fluxes  at  a  given  point  to  compute 
the  deepening  or  shallowing  of  the  mixed  layer,  without  regard 
to  horizontal  or  vertical  advective  effects.  Vertical 
advection  can  be  incorporated  into  a  one -dimensional  model,  as 
done  by  Luan  (1993)  with  PWP,  and  by  Adamec  et  al .  (1981)  and 
Muller  et  al.  (1984)  with  Garwood,  but  this  modification  was 
not  made  to  either  model  used  here  due  to  the  uncertainty 
involved  in  estimating  its  magnitude  and  variability. 
Advection,  of  course,  is  likely  to  be  of  significance  in  a 
coastal  region  such  as  this,  where  boundary  currents  and 
upwelling  effects  are  present.  These  influences,  as  well  as 
internal  forces  from  below  the  mixed  layer,  were  acknowledged 
to  be  a  likely  source  of  model  error,  but  it  was  expected  that 
they  would  be  roughly  quantifiable  through  comparison  of  the 
model  temperature  patterns  with  observations. 

The  same  forcing  was  used  to  run  each  model.  Computation 
of  these  air/sea  fluxes  is  discussed  in  Chapter  IV.  Both  use 
a  vertical  resolution  of  1  meter  and  a  time  step  of  1  hour, 
for  reasons  discussed  in  Chapter  V.  The  following  subsections 
will  present  the  pertinent  aspects  of  the  two  models, 
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demonstrating  that,  although  both  stem  from  the  integrated 
model  of  Kraus  and  Turner  (1967)  ,  they  are  fundamentally 
different  in  their  approaches  to  mixing  and  entrainment  of 
water  from  below  the  mixed  layer.  The  Garwood  model  is  based 
on  the  turbulent  kinetic  energy  equation  within  the  mixed 
layer,  while  PWP  uses  a  Richardson  number  criterion  to 
determine  when  mixing  to  another  level  should  occur.  The 
reader  is  referred  to  the  original  papers  for  details  not 
contained  in  the  following  paragraphs,  i.e.,  Garwood  (1977) 
and  Price  et  al .  (1986). 

A.  THE  GARWOOD  (1977)  MODEL 

The  basis  for  this  model  is  the  budget  for  turbulent 
kinetic  energy  (TKE)  within  the  mixed  layer.  It  has  its  roots 
in  the  basic  Kraus  &  Turner  model,  but,  like  others  developed 
since,  includes  a  parameterization  of  viscous  dissipation. 
Additionally,  it  is  unique  in  its  recognition  of  the  non¬ 
isotropic  nature  of  mixed  layer  turbulence,  and  accounts  for 
this  by  breaking  up  the  TKE  equations  into  horizontal  and 
vertical  components  (McCormick  and  Meadows,  1988) .  It  was 
presented  by  Garwood  (1977)  as  a  model  capable  of  simulating 
cyclical  steady  states  by  the  mixed  layer  from  diurnal  to 
annual  time  scales. 

Conservation  of  heat  is  generalized  to  an  equation  for 
conservation  of  buoyancy  within  the  mixed  layer,  where 
buoyancy  is  defined  as 
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(1) 


gr(Pp-g) 

Po 


and 


p=p0  [l-a  (f-T0)  +p  (s-s0)  ]  .  (2) 

The  tilde  represents  a  total  instantaneous  value  and  the 
subscript  0  indicates  an  arbitrary  representative  value.  T  is 
temperature,  s  is  salinity,  p  is  density,  and  a  and  0  are  the 
thermal  and  haline  expansion  coefficients,  respectively. 
Entrainment  is  derived  from  a  solution  of  the  total  turbulent 
kinetic  energy  equation, 


_1  d  ( u2+v2+w2)  _ 
2  dt 


r  —  dU  —  dVi  d  r  * 

-  +vw^—]  ~^bw~ -r —  [  w( 

dz  az  az 


u2  +  v 2- 


W- 
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where  upper  case  letters  denote  mean  quantities  and  lower  case 
represents  fluctuating  components.  Overbars  indicate 
instantaneous  quantities  that  have  been  averaged  over  time. 
The  letters  u,  v,  and  w  represent  eastward,  northward,  and 
upward  velocities  and  e  is  viscous  dissipation.  In  this 
equation,  the  rate  of  entrainment  depends  on  the  rate  of 
supply  of  energy  from  above  the  entrainment  zone,  which  is  the 
third  term  on  the  right  hand  side  of  equation  (3)  ,  computed  at 
depth  h.  The  entrainment  mechanism  is  theorized  to  be  of 
Benjamin's  (1963)  class  C,  in  which  local  Helmholtz 
instabilities  act  to  advect  packets  of  denser  water  up  into 
the  mixed  layer.  An  entrainment  time  scale,  re,  is  defined  as 
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the  time  needed  to  transport  turbulent  energy  to  the 
entrainment  zone.  This  is  given  by 

x  ^-a1iAw2'rin  14) 

where  ax  is  a  constant  of  proportionality,  <w2>1/2  is  the  rms 
vertical  velocity  scale,  and  the  angle  brackets  represent 
vertical  integration  over  the  entrainment  zone. 

Assuming  momentum  and  buoyancy  transport  below  the 
entrainment  zone  are  negligible,  the  mean  buoyancy  and 


momentum  equations 


dB  _  dbw  ^  ag  q 
dt  dz  pQCp 


dU  fv_  duw 

at  dz 


dV_  _  fn_  dvw 
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can  be  integrated  across  the  entrainment  zone  of  thickness  6 
to  yield  the  following  so-called  jump  conditions 


-bw(  -h)  =A B- 
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-vw(-h)  =Av||, 


18 


where  the  A' s  represent  a  jump  from  above  to  below  the 
entrainment  zone.  Garwood  postulates  that  the  critical  factor 
in  determining  the  entrainment  rate  is  the  quantity 


P= 


DW 


4- 

dz 


E 

2 


(11) 


where  E  is  the  TKE  (u^+v2+w2)  and  the  relationship  represents 
the  ratio  of  buoyancy  flux  to  convergence  of  energy  flux. 
This,  together  with  equation  (4)  produces  the  Garwood 
entrainment  equation 


P(-h) 


hJ5w{  -h)  =m 
(i^)1/2(E) 


(12) 


where  m4  becomes  the  first  of  five  empirically  derived 
dimensionless  constants  in  the  model.  With  the  introduction 
of  two  new  unknowns,  <E>  and  <w2>,  the  prediction  of  the  upper 
ocean  thermal  profile  requires  closure  through  mean  turbulent 
field  modeling  of  the  vertically  integrated  turbulent  kinetic 
energy  component  equations,  plus  the  bulk  equations  for 
buoyancy  and  momentum  (Garwood,  1977) .  Vertical  integration 
of  equations  (5) ,  (6)  ,  and  (7)  across  the  depth  of  the  mixed 
layer,  with  assumptions  of  negligible  vertical  fluxes  below 
the  mixed  layer  and  homogeneity  within  the  layer,  yields  the 
bulk  relationships  for  buoyancy  and  momentum  in  the  model. 
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.The  Heaviside  step  function  A  is  zero  in  the  case  where  the 
mixed  layer  is  shallowing,  and  1  otherwise.  Again,  horizontal 
advection  of  buoyancy  and  momentum  are  neglected. 

Garwood  defines  a  convective  turbulent  eddy  time  scale, 
tx,  proportional  to  mixed  layer  depth  divided  by  rms  turbulent 
velocity,  and  a  rotational  time  scale,  r2,  equal  to  the 
inverse  Coriolis  parameter.  These  time  scales  are  combined 
according  to 


(16) 


to  produce  the  Garwood  version  of  the  turbulent  dissipation 


rate,  or 
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where  R0  =  u*/hf  is  a  Rossby  number  for  the  mixed  layer  and  1% 
and  m5  are  two  new  empirically  derived  constants. 

The  bulk  equation  for  pressure  redistribution  R  used  in 
the  model  is 
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(18) 


R =nu(E)3/  2  -  2m~(E)1/2/ w-  > . 

c  *- 

With  an  assumption  that  the  mean  layer  velocity  is 
proportional  to  u, ,  the  surface  friction  velocity,  the  final 
term  of  the  model  may  be  generated.  This  is  the  net  wind 
generated  shear  production  rate,  or 

r  -  dV„  3  ,  1  „  2 
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The  preceding  development  results  in  the  final  equations 
which  represent  the  Garwood  (1977)  model.  These  are 
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where  C  =  U  +  iV,  mx  through  m5  are  empirical  constants,  and 
the  other  quantities  are  as  previously  defined.  Equation  (20) 
represents  the  entrainment  buoyancy  flux,  while  equations  (21) 
and  (22)  are  the  horizontal  and  vertical  components  of 
turbulent  kinetic  energy  within  the  mixed  layer,  respectively. 
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B.  THE  PRICE,  WELLER  &  PINKEL  (1986)  MODEL 

This  model  also  benefits  from  the  assumption  of  thermal, 
haline  and  velocity  uniformity  within  the  mixed  layer.  Its 
fundamental  difference,  following  from  Pollard  et  al .  (1973), 
is  the  way  in  which  it  determines  when  mixing  to  another 
vertical  level  should  occur. 

The  PWP  model  begins  with  the  one -dimensional  heat  and 
momentum  equations,  as  follows 


ar_  1  dQ 

dc~  P0 cp  dz 


(23) 


dC  p0  dz 


(24) 


£'—fu--±-Zp,  (25) 

dt  p0  dz 

where  the  Gx  and  Gy  represent  the  components  of  the  shear 
stress  in  the  water  column,  which  at  the  surface  equals  the 
wind  stress.  Equations  (24)  and  (25)  are  modified  in  the 
version  of  the  PWP  model  used  here  to  include  a  frictional 
term,  as  done  in  Rosenfeld  (1988) .  The  value  assigned  to  the 
friction  parameter  is  discussed  in  Chapter  V. 

As  in  Garwood  (1977) ,  density  is  calculated  using  a 
linearized  equation  of  state,  similar  to  equation  (2)  . 
Absorption  of  solar  insolation  is  handled  in  much  the  same  way 
as  it  is  in  the  Garwood  model.  Both  absorb  about  half  of  the 
incoming  irradiance  in  the  upper  meter  of  water,  with  an 
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exponential  decay  below  that.  PWP  additionally  breaks  the 
solar  energy  into  long  and  shortwave  components  and  absorbs 
them  according  to  different  vertical  decay  scales.  A  brief 
discussion  of  water  type  assumptions  made  to  determine  the 
absorption  coefficients  is  contained  in  the  model 
sensitivities  section  of  Chapter  V. 

Mixing  is  modeled  as  in  the  dynamic  instability  model  of 
Price  et  al .  (1978).  It  takes  place  in  such  a  way  as  to 
satisfy  conditions  for  static,  mixed  layer  and  shear  flow 
stability,  which  are  represented  by  equations  (26) ,  (27) ,  and 
(28) ,  respectively. 

i£*0  (26) 

dz 
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uppermost  meter  of  the  water  column  at  each  time  step. 
Densities  are  then  computed  and  mixed  to  achieve  static 
stability.  Next,  wind  stress  is  absorbed  and  a  value  of  Rb  is 
computed  and  compared  to  the  critical  value  of  0.65. 
Entrainment  of  denser  water  from  below  the  mixed  layer  takes 
place  as  necessary  to  achieve  condition  (26).  Price  et  al . 
(1986)  also  include  the  shear  flow  stability  requirement  to 
account  for  the  fact  that  their  observations  indicate  a  smooth 
transition  layer  below  the  mixed  layer,  as  opposed  to  the 
sharp  jump  assumed  by  other  models,  and  because  shear 
instability  is  likely  at  levels  of  strong  stratification.  Rg 
is  calculated  only  within  this  transition  layer  and  mixing 
takes  place  until  it  is  above  the  critical  value  of  0.25 
throughout  this  region.  This  has  the  effect  of  smoothing  out 
the  transition  from  mixed  layer  to  the  region  below. 


24 


IV.  DATA 


A.  SOURCES 

1.  OASIS  MOORING 

A  modified  ATLAS  (Automated  Temperature  Line 
Acquisition  System)  mooring  (Milburn  and  McClain,  1986) , 
located  at  36°45'N,  122°01'W  (see  Figure  1) ,  was  the  primary 
source  of  meteorological  and  oceanographic  data  used  in  this 
study.  The  modification,  done  at  MBARI,  includes  a  controller 
known  as  OASIS,  Ocean  Acquisition  System  for  Interdisciplinary 
Science  (Chavez  et  al.,  1991) .  The  mooring,  designated  as  Ml, 
is  owned  and  maintained  by  MBARI,  with  Dr.  Francisco  Chavez 
acting  as  supervising  scientist.  It  is  schematically  depicted 
in  Figure  3a,  with  the  details  of  the  buoy  itself  shown  in 
Figure  3b.  Above  the  waterline,  it  is  equipped  with  sensors 
for  air  temperature,  humidity,  wind  speed  and  direction  and 
photosynthetically  active  radiation  (PAR) .  Below  the  surface 
is  a  thermistor  chain  with  sensors  at  10,  20,  40,  60,  80,  100, 
150,  200,  250,  and  300  meters.  CTD  sensors  are  located  at  the 
surface,  10,  and  20  meters.  Other  sensors  include  a 
fluorometer  and  a  transmissometer ,  although  these  were  not 
specifically  used  in  this  study.  All  of  these  data  are 
collected  at  a  10  minute  interval  and  transmitted  to  MBARI  via 
packet  radio  telemetry.  Winds  are  measured  3.8  meters  above 
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the  sea  surface  at  a  frequency  of  2  Hz  and  averaged  over  one 
minute.  The  mooring  is  also  equipped  with  an  acoustic  Doppler 
current  profiler  (ADCP) ,  which  samples  at  a  15  minute 
interval.  Data  quality  for  all  parameters  of  interest  to  this 
study  was  excellent  throughout  the  period  with  no  gaps  or 
systematic  errors,  with  the  exception  of  the  sea  surface 
temperature  sensor  on  the  ATLAS  mooring,  which  read 
consistently  high  by  a  few  degrees.  Sea  surface  temperatures 
used  were  therefore  obtained  from  the  CTD  sensor  mounted  on 
the  mooring  at  the  surface.  The  time  response  of  this  sensor 
is  on  the  order  of  seconds,  whereas  the  time  response  for  the 
subsurface  sensors  is  on  the  order  of  minutes . 

A  possible  error  in  the  data  was  in  the  measurement  of 
wind  direction.  This  is  determined  on  board  the  buoy  by 
combining  a  compass  measurement  of  the  buoy  orientation  with 
an  anemometer  vane  measurement  of  the  wind  direction  relative 
to  the  buoy.  Later  in  the  mooring  deployment,  the  compass 
failed.  While  examination  of  the  data  used  here  show  no 
conclusive  evidence  of  a  problem  with  wind  direction, 
comparisons  with  data  from  other  nearby  land  stations  suggest 
that  wind  direction  could  be  off  by  up  to  30°.  Since  it  is 
only  the  strength  and  rotation  rate  of  the  wind  that  is 
important  in  a  one -dimensional  analysis,  this  possible  error 
should  not  present  a  problem  for  this  study.  Specifications 
and  details  of  operation  for  the  OASIS  mooring  are  contained 
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in  the  Pacific  Marine  Environmental  Laboratory  reference 
manual  (1992) . 

2.  NPS  Profiler  Site 

Shortwave  and  longwave  irradiance  data  used  in  this 
study  were  collected  at  the  Naval  Postgraduate  School's  (NPS) 
meteorological  station  located  on  Fort  Ord  as  shown  in  Figure 
1.  This  site,  located  23  kilometers  east- southeast  of  the 
mooring,  is  equipped  to  measure  a  variety  of  meteorological 
parameters  at  2  minute  resolution.  However,  only  the 
irradiance  data  was  of  practical  use  due  to  distance  from  Ml 
and  availability  of  other  necessary  data  from  the  mooring 
itself.  Shortwave  irradiance  is  measured  with  an  Eppley 
Precision  Pyranometer  (Model  PSP)  and  longwave  irradiance  is 
measured  with  an  Eppley  Precision  Infrared  Radiometer  (Model 
PIR)  .  Specifications  for  both  are  contained  in  the  Eppley 
Laboratory  reference  manual  (Eppley  Laboratory,  1971) . 

B.  DESCRIPTION  OF  ATMOSPHERIC  AND  OCEANIC  DATA  AND  DERIVED 

QUANTITIES 

1.  Synoptic  Weather  Description 

During  most  of  the  period  of  study,  Monterey  Bay  was 
under  the  influence  of  the  typical  summer  high  pressure  system 
over  the  eastern  North  Pacific  Ocean.  A  1030  mb  high  centered 
near  45 °N,  160 °W  can  be  seen  in  Figure  4,  which  indicates  the 
surface  pressure  pattern  at  1200Z  on  3  September,  and  which  is 
representative  of  the  period  l  to  8  September.  Generally 
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clear  skies  inland  with  some  low-level  stratus  offshore  can  be 


seen  in  the  imagery  from  this  time  period.  The  result  of  this 
pattern,  which  persists  for  long  intervals  in  this  region,  is 
light  and  variable  synoptic  scale  winds,  with  subsidence 
associated  with  the  eastern  portion  of  this  system  creating 
the  usual  low-level  marine  inversion  over  the  cool  upwelled 
waters  of  the  coastal  region.  This  marine  layer  is  the  source 
of  the  frequent  night  and  morning  fog,  which  occurs  to  varying 
degree  on  each  of  the  days  of  the  study  period. 

The  clear  inland  weather  associated  with  this  fair 
weather  system,  combined  with  the  weak  large  scale  forcing, 
allowed  for  extensive  heating  in  the  Salinas  Valley.  High 
temperatures  in  King  City,  near  the  head  of  the  valley, 
reached  to  near  100 °C  during  most  of  the  period.  This  heat, 
in  the  presence  of  the  light  synoptic  scale  winds,  allowed  for 
near  classical  development  of  sea  breeze  conditions  in  the 
Monterey  Bay. 

During  the  9th  and  10th  of  September,  the  high  pressure 
weakened  to  1023  mb  and  gave  way  along  the  coast  to  low 
pressure  centered  off  northern  California.  This  produced 
significant  cloudiness  in  the  study  area.  The  surface 
pressures  from  10  September  are  shown  in  Figure  5.  Synoptic 
scale  flow  remained  weak  during  this  time,  but  the  cloudiness 
greatly  reduced  the  heating  onshore  and  thus  the  sea  breeze 
signature,  as  will  be  seen  in  the  following  sections. 
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2 .  Hourly  Data 

The  following  subsections  present  the  raw  wind,  sea 
surface  temperature,  air  temperature,  humidity,  and  10,  20  and 
40  meter  sea  temperatures  collected  at  mooring  Ml.  Radiation 
data  will  be  presented  in  the  next  section,  along  with  a 
discussion  of  its  use.  Corresponding  figures  are  located  at 
the  end  of  this  chapter.  The  data  are  shown  with  the  full 
temporal  resolution  of  10  minutes.  However,  it  is  noted  here 
that  hourly  averages  of  wind  components,  surface  and  air 
temperatures,  and  humidity  were  computed  and  used  to  generate 
the  model  forcing  for  reasons  discussed  in  Chapter  V.  A 
comparison  of  raw  and  hourly  averaged  data  indicated  that  all 
of  the  key  features  were  maintained  in  this  process, 
particularly  the  timing  and  duration  of  each  diurnal  event. 

a.  Winds 

U  and  V  wind  components  are  shown  in  Figure  6a, 
while  total  wind  speed  is  shown  in  Figure  6b.  The  east /west 
component  dominates  the  signal  on  most  days.  Readily  apparent 
is  the  strongly  diurnal  variation  in  the  magnitude  of  the 
winds.  Wind  magnitude  tends  to  be  quite  light  at  night  and 
into  early  morning,  averaging  about  2  m/s  over  the  11  day 
period.  At  an  average  time  of  1737Z  (1037  PDT)  ,  the  sea 
breeze  front  reaches  the  mooring  location  and  winds  accelerate 
rapidly,  reaching  a  peak  of  6  to  12  m/s  by  the  mean  time  of 
0029Z  (1729  PDT)  ,  averaged  over  all  11  days  of  the  period. 
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The  peak  wind  of  the  period  occurred  at  0116Z  on  the  6th, 
which  corresponds  to  1816  PDT  on  the  5th,  when  the  wind 
exceeded  12  m/s.  A  plot  of  true  wind  direction  (Figure  6c) 
shows  that  the  rapid  increase  in  wind  speed  associated  with 
the  onset  of  the  sea  breeze  corresponds  to  the  abrupt  shifting 
of  the  direction  toward  the  southeast,  or  about  130  degrees 
true,  as  would  be  expected  by  the  orientation  of  the  Salinas 
Valley.  In  general,  the  winds  taper  off  more  gradually  in  the 
late  afternoon  and  are  replaced  with  the  light,  variable 
winds,  with  little  or  no  tendency  to  favor  an  offshore  land 
breeze  toward  the  west.  This  offshore  wind  is  better 
developed  at  the  profiler  site.  Another  view  of  the  winds  is 
shown  in  Figure  6d,  a  feather  plot  of  hourly  averaged  wind 
vectors,  which  shows  that  the  wind  is  generally  blowing  toward 
the  southeast,  but  intensifies  most  afternoons  from  a  mean 
direction  of  310°T. 

It  should  be  noted  that  the  strongest  sea  breeze 
days  were  6,  7,  and  8  September,  when  the  inland  and  coastal 
weather  was  the  clearest  and  maximum  heating  occurred  in  the 
valley.  Correspondingly,  the  wind  magnitudes  observed  at  the 
mooring  were  extremely  light  on  9  through  11  September,  when 
extensive  cloudiness  was  present  both  in  the  bay  and  at  the 
profiler  site,  as  well  as  further  inland.  The  winds  on  these 
days  only  reached  maxima  of  5  to  6  m/s. 

The  winds  can  clearly  be  characterized  as  being 
overwhelmingly  dominated  by  variations  in  the  diurnal 
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frequency  band,  although  to  a  reduced  degree  during  this 
latter  period.  Therefore,  the  momentum  fluxes  computed  from 
these  winds  will  be  strongly  diurnal.  The  phase  relationship 
of  this  diurnal  wind  stress  pattern  to  the  diurnal  heating 
will  be  discussed  later  in  this  chapter. 

Jb.  Air  Temperature 

As  would  be  true  in  most  areas,  air  temperatures 
also  exhibit  a  strongly  diurnal  variation,  although  the 
magnitude  of  the  day/night  swing  is  small  compared  to  the 
nearby  land  data.  Temperatures  observed  at  Ml  in  early 
September  are  seen  in  Figure  7.  The  mean  air  temperature  over 
the  period  was  14.41°C,  with  a  standard  deviation  of  0.83. 
The  average  high  temperature  was  15.22°C,  and  the  low  averaged 
13 . 50°C.  An  average  diurnal  temperature  swing  of  l.7l°C 
occurs  over  the  period,  although  of  note  is  the  3.25°C  range 
seen  from  maximum  to  minimum  on  the  6th.  Mean  time  of  warmest 
temperatures  for  the  day  is  0115Z  (1815  PDT) ,  with  a  range  of 
2226Z  to  0353Z.  Daily  lows  occur  between  1206Z  and  1928Z, 
with  an  average  of  1658Z  (0958  PDT)  . 

c.  Humidity/Dew  Point 

As  mentioned  previously,  the  ATLAS  mooring 
measures  relative  humidity.  These  data  are  shown  in  Figure  8. 
For  the  purposes  of  computing  fluxes,  as  will  be  discussed  in 
the  next  section,  a  conversion  to  dew  point  temperatures  was 
made.  These  are  plotted  along  with  the  air  temperatures  in 
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Figure  7.  As  would  be  expected,  the  air  was  quite  moist 
throughout  the  period,  with  the  exception  of  the  late 
afternoon  of  the  5th  (about  0000Z  on  the  6th)  ,  when  the 
relative  humidity  dipped  briefly  to  76%  with  a  strong, 
relatively  dry  sea  breeze  due  to  the  presence  of  drier  air 
offshore.  Otherwise,  the  humidity  is  consistently  above  90% 
with  the  air  at  or  near  saturation  on  several  evenings  and 
mornings,  as  well  as  most  of  the  9th  through  the  llth,  related 
to  fog  and/or  low  clouds. 

d.  Sea  Surface  Temperature 

Ocean  surface  temperatures  as  measured  by  the 
surface  CTD  on  the  OASIS  mooring  are  presented  in  Figure  9 . 
The  diurnal  cycle  is  again  the  dominant  feature  of  this 
pattern.  Over  the  11  day  period,  the  sea  surface  temperature 
averages  14.64°C,  with  a  standard  deviation  of  0.51.  Daily 
maxima  occur  between  2027Z  and  0223Z,  with  an  average  time  of 
2303Z  (1603  PDT)  ,  and  have  a  mean  value  of  15.30°C.  The 
minimum  temperature  averages  14.26°C  and  occurs  at  1334Z  (0634 
PDT)  in  the  mean.  Thus,  the  surface  temperature  minima  and 
maxima  lead  those  of  the  air  temperature  by  3:24  and  2:12  h, 
respectively.  The  magnitude  of  the  diurnal  variation  in  sea 
surface  temperature  data  will  be  of  importance  in  later  model 
result  discussions.  This  day  to  night  swing  ranges  from 
0.57°C  to  1.55°C  and  averages  0.96°C. 
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e.  Temperatures  at  10,  20,  and  40  Meters 

Ocean  temperatures  from  the  surface,  10,  20  and  40 
meter  thermistors  are  shown  in  Figure  10.  At  10  meters,  the 
diurnal  pattern  is  still  quite  apparent  and  dominates  the 
temperature  signal.  Occasions  on  which  the  mixed  layer 
deepens  to  or  below  this  depth  are  seen  here  as  convergences 
of  the  surface  and  10  meter  traces.  It  is  clear  that  this 
occurs  on  several  nighttime  cooling  cycles,  but  not  at  all 
during  the  latter  portion  of  the  period  when  winds  were  light. 
Direct  measurements  of  mixed  layer  depth  were  not  available 
during  the  study  period  {with  the  exception  of  a  series  of  CTD 
casts  made  on  9  September  on  a  MBARI  cruise) .  The  20  meter 
temperatures  show  a  signal  dominated  by  much  higher  frequency 
variation  than  diurnal,  driven  by  internal  waves  rather  than 
surface  forcing.  During  September  4th,  it  appears  that  the 
mixed  layer  actually  deepened  briefly  to  this  depth.  At  40 
meters,  a  similarly  high  frequency  temperature  pattern  is 
seen,  although  quite  damped  relative  to  20m.  Of  note  is  a 
slight  downward  trend  in  the  10  and  20  meter  data  over  the 
period.  This  will  be  discussed  more  quantitatively  with 
respect  to  advective  effects  in  Chapter  VI. 

C.  MIXED  LAYER  CHARACTERISTICS  AT  Ml 

The  data  described  above  depict  an  ocean  mixed  layer  that, 
like  in  other  areas,  undergoes  a  strongly  diurnal  variation  in 
thermal  structure.  The  mostly  clear  weather  permits 
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significant  heating  and  resultant  stratification  during  the 
day,  with  relatively  calm  conditions  and  convective  cooling  at 
night.  The  mixed  layer  depth  is  generally  less  than  20  meters 
and  often  shoals  to  less  than  10  meters,  particularly  when  the 
winds  are  light  and  the  insolation  is  strong.  The  wind 
pattern  described  is  also  varying  diurnally,  with  near  calm 
conditions  being  rapidly  replaced  with  winds  of  8  to  12  m/s. 
The  peaks  of  wind  stress  occur  in  the  late  afternoon  to  early 
evening  each  day,  approximately  coincident  with  the  peak  in 
air  temperature  but  well  after  that  of  solar  insolation. 
Thus,  the  afternoons  are  characterized  by  the  strongest 
turbulence  generation  due  to  the  wind,  but  also  fairly  strong 
buoyant  damping  from  the  continuing  downward  heat  flux.  The 
winds  die  down  fairly  rapidly  as  darkness  sets  in,  reducing 
the  shear  generation  of  turbulence,  but  at  the  same  time, 
convective  cooling  is  creating  buoyant  turbulence.  By  early 
morning,  buoyant  damping  caused  by  the  rising  sun  can  quickly 
shallow  the  mixed  layer,  since  the  winds  are  typically  at 
their  weakest  at  that  time.  The  next  chapter  will  present  a 
discussion  of  model  sensitivities,  followed  by  the  model 
results  themselves.  A  discussion  of  their  ability  to  handle 
mixed  layer  behavior  in  this  relatively  complex  environment 
will  follow  in  Chapter  VI. 
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D .  SURFACE  FLUXES 


1.  Cloud  Cover  Estimation 

No  direct  measurements  of  cloud  cover  were  made  during 
the  study  period.  In  order  to  provide  reasonable  cloud  cover 
information  for  the  purpose  of  generating  valid  surface 
fluxes,  it  was  necessary  to  estimate  cloud  cover  from 
available  indirect  methods.  To  achieve  this,  insolation  and 
longwave  downward  irradiance  data  from  the  nearby  NPS  profiler 
site  and  photosynthetically  active  radiation  (PAR)  data  from 
the  mooring  were  integrated  to  create  the  best  possible 
estimation  of  percent  cloud  cover. 

These  data  are  shown  in  Figures  lla  and  b.  Insolation 
data  was  compared  with  that  observed  on  7  September  (an 
extremely  clear  day)  and  10  September  (a  completely  cloud 
covered  day)  to  obtain  an  approximation  of  cloud  cover  in 
eighths  for  the  daytime  periods.  The  presence  of  fog  or  low 
clouds  is  evident  in  the  data  during  the  intervals  when  the 
value  of  longwave  irradiance  jumps  to  about  400  W/m2,  such  as 
1000  through  1900  GMT  on  the  1st,  and  remains  fairly  steady  at 
that  level.  This  was  very  useful  for  confirming  and  adjusting 
the  cloud  cover  during  the  daylight  hours,  but  also  served  to 
provide  a  means  for  rather  accurate  assessment  of  clouds  at 
night.  A  steady  400  W/m2  is  consistently  indicative  of 
complete  low-level  stratus  or  fog.  These  downward  irradiance 
data  demonstrate  that  the  cloud  cover  is  frequently  either 
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100%  or  0%,  as  there  are  only  relatively  brief  occasions  when 
an  intermediate  value  is  measured. 

The  combination  of  the  downward  irradiance  and 
incoming  shortwave  radiation  produced  an  estimate  of  cloud 
cover  at  the  profiler  site.  The  PAR  data  were  then  used  to 
assess  the  validity  of  this  cloud  pattern  at  the  mooring.  In 
general,  the  PAR  sensor  confirmed  the  mainly  sunny  pattern, 
with  some  clouds  at  the  end  of  the  period,  seen  at  the 
profiler.  On  the  9th  (peak  centered  prior  to  10/00Z  in  Figure 
lib) ,  however,  it  was  clear  chat  a  much  cloudier  day  occurred 
at  Ml  than  at  the  profiler.  This  information  was  then  used  to 
increase  the  estimate  of  cloud  cover  at  the  buoy  on  that  day. 

2.  Model  Forcing 

As  mentioned  previously,  the  same  program  was  used  to 
generate  the  heat  and  momentum  fluxes  for  both  models.  Hourly 
averages  of  U  and  V  wind  components,  air  temperature  and  dew 
point,  sea  surface  temperature,  and  cloud  cover  were  computed 
from  the  data,  and  used  as  the  input  for  the  flux  generation. 
Latent  and  sensible  heat  fluxes,  as  well  as  wind  stresses, 
were  computed  in  accordance  with  the  methods  of  Large  and  Pond 
(1981) .  Net  longwave  irradiance  and  incoming  shortwave 
irradiance  were  computed  using  the  formulations  of  Husby  and 
Seckel  (1975) .  These  calculations  were  combined  into  a 
forcing  program  by  R.W.  Garwood  of  NPS.  Adjustments  were  made 
for  latitude  and  longitude  and  the  code  was  altered  to  allow 
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hourly,  vice  three-hourly,  inputs  of  the  meteorological  values 
in  order  to  better  reflect  the  diurnal  time  scale  mixed  layer 
evolution.  Sensible,  latent,  and  infrared  heat  fluxes  were 
combined  into  a  net  heat  loss  value.  This,  together  with  the 
solar  insolation  values  and  wind  stresses,  was  passed  to  the 
models  for  forcing,  with  adjustments  made  for  different  format 
and  unit  requirements  in  each  model.  Plots  of  all  the  forcing 
values  computed  and  passed  to  the  models  are  contained  in 
Figures  12a,  b,  c,  and  d. 

The  first  of  these  is  computed  wind  stress,  using  the 
relationship 


T  =  PaCDU.2*  (29) 

where  CD  is  the  drag  coefficient  from  Large  and  Pond  (1981) 
and  pa  is  the  density  of  air  (1.23  kg/m3) .  Figure  12b  shows 
the  solar  insolation  pattern,  adjusted  from  the  ideal, 
cloudless  sine  wave  pattern  to  reflect  the  cloud  cover 
influences  estimated  as  described  in  the  last  subsection. 
Figure  12c  presents  the  individual  heat  loss  terms,  with 
positive  values  indicating  heat  lost  from  the  ocean  to  the 
atmosphere,  along  with  the  total  hourly  heat  loss  passed  to 
the  models  for  forcing.  From  this  last  figure,  it  is  clear 
that  the  largest  term  in  the  net  loss  is  the  infrared 
radiation,  Qb.  This  reaches  a  minimum  of  about  3  0  W/m2  during 
cloudy  periods  and  a  maximum  of  about  75  W/m2  on  clear  days. 
The  sensible  heat  flux  can  be  seen  to  generally  be  near  0, 
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reaching  a  maximum  loss  of  7  W/m2  late  on  the  7th  and  a 
minimum  of  -25  W/m2  early  on  the  sixth  (negative  values 
indicate  a  gain  of  heat  by  the  ocean) ,  when  a  brief  surge  of 
relatively  dry,  warm  air  passed  through  the  mooring  area.  The 
latent  heat  flux  is  obviously  closely  tied  to  the  wind  stress 
pattern  shown  in  Figure  12a.  Finally,  a  summary  plot  (Figure 
I2d)  shows  the  net  heat  flux  computed  by  the  forcing  routine, 
demonstrating  the  dominance  of  the  short  wave  solar  heating. 

A  statistical  summary  of  the  wind  stress  and  net  heat 
flux  is  offered  here  for  later  use  in  the  discussion  of  the 
model  results.  Wind  stress  averaged  0.340  dynes/cm2  over  the 
period,  with  a  standard  deviation  of  0.407.  The  average 
minimum  of  0.044  occurred  at  173 0Z  (1030  PDT)  in  the  mean. 
The  maximum  averaged  0.942  dynes/cm2  and  occurred  typically  at 
0029Z  (1729  PDT) ,  although  on  the  stronger  wind  days  in  the 
middle  of  the  period,  this  peak  occurred  closer  to  1840  PDT. 
The  maximum  downward  heat  flux  each  day  had  a  mean  value  of 
650  Watts/m2.  The  timing  of  this  peak  was  very  regular  at 
2011Z  (1311  PDT) . 

Black  body  irradiance  from  the  sea  surface  was 
combined  with  measured  longwave  down  data  as  a  check  on  the 
validity  of  the  computed  net  infrared  heat  loss,  the  largest 
term  in  the  total  heat  loss.  Values  were  in  good  agreement 
(within  a  maximum  error  of  20  W/m2)  .  Varying  heat  loss  within 
the  range  of  reasonable  values  showed  the  insensitivity  of  the 
models  to  this  parameter,  relative  to  wind  and  insolation. 
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Configuration:  Courtesy  of  Francisco  Chavez,  MBARI. 
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Figure  3b.  Schematic  Diagram  of  the  ATLAS  Buoy:  Courtesy 
of  Francisco  Chavez,  MBARI. 
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Figure  4.  Surface  Pressure  Chart  From  4  September,  1992. 
Monterey  Bay  is  Designated  With  a  *.  Solid  Lines  are 
Surface  Pressures  in  Millibars  and  Dashed  Lines  Indicate 
1000  to  500  mb  Thickness. 


Figure  5.  Surface  Pressure  Chart  From  10  September,  1992. 
Same  Line  Convention  as  in  Figure  4. 


Figure  6a.  Full  Resolution  U  and  V  wind  Components 
Measured  at  Ml.  The  V  Component  is  Dashed.  This  and 
All  Subsequent  Data  and  Model  Plots  are  in  Greenwich  Mean 
Time.  Local  Time  May  be  Obtained  by  Subtracting  7  Hours. 
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September,  1992 
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Figure  6c.  True  Wind  Direction  at  Ml.  Values  Indicate 
the  Direction  the  Wind  Vector  was  Pointing.  Periods  of 
Large  Variation  Occur  During  Light  Winds. 
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September,  1992 


o  S00j60a 

Figure  7.  Air  Temperature  (Solid)  and  Dew  Point  (Dashed) 
at  Ml.  Dew  Points  Were  Converted  From  Relative  Humidities 
Measured  by  the  Buoy. 
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16.5 


September,  1992 


Figure  10.  Ocean  Temperatures  at  the  Surface,  10,  20,  and 
40  Meters  Beneath  Ml. 
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September,  1992 
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Figure  11a.  Shortwave  and  Longwave  Irradiance  Recorded  at 
the  NPS  Profiler  Site.  The  Day  of  September  is  Designated 
on  the  Time  Scale  of  Each  Plot. 
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Figure  11a.  (Continued) 
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Figure  11a.  (Continued) 
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Figure  11a.  (Continued) 
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2000 


Figure  lib.  Photosynthetically  Active  Radiation  at  Ml. 
Note  Degree  of  Cloudiness  on  the  9th. 
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Figure  12a.  Total  Wind  Stress  Computed  in  Accordance  With 
Large  and  Fond  (1981) . 
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September,  1992 


CM 


Figure  12c.  Computed  V  at  Loss  Terms.  Sensible  (Dotted), 
Latent (Lower  Solid),  Net  Infrared  (Dashed),  and  Net  Heat 
Loss  (Upper  Solid  -  the  Other  3  Summed)  are  Plotted,  with 
Positive  Values  Indicating  Heat  Lost  From  the  Ocean. 
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September,  1992 


Figure  12d.  Computed  Net  Surface  Heat  Flux.  Sign 
Convention  as  in  Figure  12c.  This  Represents  the  Solar 
Insolation  Minus  the  Net  Heat  Loss  From  Figure  12c. 
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September,  1992 


V.  MODEL  RESULTS 


A.  MODEL  SENSITIVITIES 

1.  Initial  Conditions 

In  order  to  run  the  two  mixed  layer  models,  an  initial 
density  profile  was  needed.  A  CTD  or  bathythermograph  cast 
was  not  available  for  the  beginning  of  the  1st  of  September, 
when  the  data  collection  began,  but  a  MBARI  cruise  made  on  9 
September  conducted  nine  CTD  casts.  All  of  the  casts  were 
made  within  the  bay.  The  temperatures  from  the  upper  50 
meters  of  each  cast  are  shown  in  Figure  13a.  Initial  runs  of 
the  models  were  made  using  the  cast  from  this  cruise  which 
most  nearly  matched  the  surface  temperature  with  that  of  the 
ATLAS  buoy  at  0000Z  on  1  September  1992  (designated  #1  in 
Figure  13a) .  Experimentation  with  a  variety  of  other  initial 
temperature  profiles  indicated  that  the  early  model  heat 
content  in  the  mixed  layer  was  strongly  dependent  on  the 
starting  temperature  profile.  Since  the  initial  profile  in 
the  early  runs  was  nearly  matched  with  the  data  in  surface 
temperature  and  time  of  day,  but  was  not  located  near  the 
mooring  (about  15.5  kilometers  to  the  northeast),  it  was 
decided  that  a  profile  taken  at  the  mooring  site  would  be  more 
representative  of  conditions  at  that  site.  Subsequent  model 
runs  were  started  at  2000Z  on  the  1st,  since  the  CTD  cast 
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located  at  the  buoy  (designated  #2  in  Figure  13a)  was  made  at 
that  time,  although  on  a  different  date.  A  0.325°C  mismatch 
in  surface  temperature  between  the  CTD  cast  from  the  9th  and 
the  ATLAS  data  on  the  1st  was  corrected  by  adding  that  amount 
to  the  CTD  cast  along  the  51  meters  of  that  cast  used  as  an 
initial  profile.  In  this  way,  the  stratification  of  the  water 
column  was  preserved,  and  the  most  representative  initial 
temperature  profile  possible  was  attained.  Final  runs 
presented  in  the  next  section  were  initialized  using  this 
profile,  which  is  shown  separately  in  Figure  13b. 

Experimentation  with  varying  initial  salinity  profiles 
was  also  conducted.  As  designed,  the  Garwood  model  assumes  a 
constant  salinity  with  depth,  while  PWP  allows  the  input  of 
observed  salinities.  For  comparison,  a  test  PWP  run  was  made, 
assuming  a  constant  salinity.  The  difference  in  mixed  layer 
depths  and  temperatures  produced  was  indistinguishable. 
Therefore,  no  modification  of  the  Garwood  model  to  allow 
observed  salinities  to  be  used  in  the  calculation  of  an 
initial  density  profile  was  deemed  necessary. 

2.  Friction  (PWP) 

As  mentioned  in  the  model  descriptions,  the  PWP  model 
contains  a  parameterization  of  friction,  which  is  treated  in 
terms  of  a  relaxation  time  of  the  surface  currents.  Initial 
runs  of  this  model  left  this  variable  at  its  default  value  of 
9999,  which  corresponds  to  no  friction.  It  was  decided  that 
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a  value  of  5  days  would  be  more  correct.  This  value  is  in 
agreement  with  that  used  by  Rosenfeld  (1988)  in  the  CODE 
region  of  northern  California.  In  terms  of  the  mixing  and 
thermal  structure  of  the  upper  ocean,  this  adjustment  made 
little  difference  over  the  11  day  study  interval,  as  can  be 
seen  in  Figure  14,  but  was  considered  to  be  more  accurate  than 
using  no  frictional  term  at  all. 

3 .  Time  Step 

Another  important  model  variable  to  be  established  was 
the  degree  to  which  it  was  reasonable  to  reduce  the  model  time 
step  in  the  interest  of  representing  the  short  temporal  scale 
mixed  layer  response  to  the  sea  breeze  environment.  It  was 
considered  that  running  the  models  at  a  10  minute  time  step  to 
correspond  with  the  available  data  might  be  ideal.  However, 
the  Garwood  model  presumes  that  TKE  reaches  a  steady  state 
within  the  mixed  layer  at  each  time  step.  To  test  the 
validity  of  this  assumption  for  short  time  scales,  a  second 
model  written  by  R.W.  Garwood  of  NPS,  in  which  the  complete 
unsteady  physics  of  the  turbulent  kinetic  energy  equations  is 
preserved,  was  used  to  test  the  time  needed  for  TKE  to  reach 
steady  state  after  a  step  wind  function  is  applied  to  the 
problem.  These  equations  are: 
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at  p 

where  h  is  mixed  layer  depth,  T  is  mixed  layer  temperature,  we 
is  entrainment  velocity,  E  is  total  TKE ,  W  is  the  mean 
vertical  velocity  at  the  base  of  the  mixed  layer,  u*  is  the 
surface  friction  velocity,  AT  is  the  temperature  change  at  the 
base  of  the  mixed  layer,  Q0/(pcp)  is  the  net  surface  heat 
flux,  R  is  a  pressure  redistribution  term,  D  represents 
viscous  dissipation,  f  is  the  Coriolis  parameter,  and  r  is 
wind  stress.  Figure  15  shows  the  wind  function,  an  idealized 
sea  breeze,  that  was  applied  to  this  model  and  the  resulting 
total  TKE  versus  time.  This  step  function  is  roughly 
equivalent  to  an  increase  of  wind  from  8  to  12  m/s.  It  can  be 
seen  that  an  hourly  or  greater  time  step  is  valid,  since  a 
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steady  state  can  be  reached  within  that  time.  However,  a 
shorter  interval  between  turbulent  calculations,  such  as  10 
minutes,  does  not  allow  this  steady  state  to  occur  and  thus 
prevents  proper  distribution  of  the  energy.  Therefore,  it  was 
determined  that  hourly  fluxes  should  be  computed  from  the 
hourly  averaged  ATLAS  data.  The  PWP  model  was  also  run  at  a 
one  hour  time  step  for  the  purpose  of  uniformity  between  the 
two  models,  although  it  is  noted  here  that  it  could  be  run 
with  a  shorter  time  step,  as  in  Price  et  al .  (1986)  ,  since  the 
mixing  is  not  based  on  a  turbulence  steady  state.  Since  the 
hourly  averaged  data  reflect  the  diurnal  cycle  well,  no  model 
runs  at  10  minute  time  steps  were  deemed  necessary. 

4 .  Cloudiness 

Another  factor  affecting  the  model  results  was  the 
cloud  cover.  As  discussed  previously,  a  best  estimate  of  this 
input  value  was  made  from  the  available  irradiance  data. 
Initial  runs  of  the  models,  however,  were  made  asstiming  no 
cloud  cover.  Figures  16a  and  b  show  the  Garwood  and  PWP  model 
mixed  layer  temperatures  with  and  without  cloud  cover.  It  is 
apparent  that  this  was  not  an  especially  important  factor 
through  most  of  the  study  period  when  the  weather  was 
primarily  clear,  except  for  periods  of  fog  at  night.  However, 
there  is  a  significant  reduction  in  the  warming  trend  of  the 
models  in  the  later  period  during  which  extensive  cloudiness 
was  present. 
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5.  Absorption  Coefficients 

Two  final  tuneable  parameters  in  the  PWP  model  are  the 
values  for  longwave  and  shortwave  absorption  coefficients. 
The  net  incoming  solar  radiation  is  absorbed  in  the  ocean 
according  to 


I(z)  =1(0)  (Ij^e  X;-Ise  is] 

where  1  and  s  designate  longwave  and  shortwave  components  of 
insolation  (Rosenfeld,  1988).  Values  of  =  .6  and  Is  =  .4 
were  used  as  in  Price  et  al.  (1986)  .  Coastal  type  III 
(Jerlov,  1976)  was  assumed  as  in  the  CODE  area,  with  values  of 
=  1.4  m'1  and  X2  =  7.9  m*1. 

Garwood  absorption  is  treated  similarly,  with 
exponential  decay  of  absorption  with  depth  below  l  meter,  but 
the  incoming  radiation  is  not  broken  into  long  and  shortwave 
components.  Figure  17  gives  the  absorption  profiles  with 
depth  for  the  two  models. 

The  following  results  describe  the  model  performance 
with  all  of  the  preceding  adjustable  parameters  set  as 
described.  They  represent  our  best  estimate  of  the  actual 
absorption,  cloudiness,  friction,  and  initial  conditions  as 
discussed  above  and  computed  at  a  one  hour  time  step  at  1 
meter  vertical  resolution. 
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B.  FINAL  MODEL  RUNS 


1 .  Garwood 

Mixed  layer  temperatures  and  depths  computed  by  the 
Garwood  model,  initialized  at  2000Z  on  the  1st  with  the 
temperature  profile  shown  in  Figure  13b,  (located  at  the 
mooring) ,  are  presented  in  Figures  18a  and  b.  Also  in  Figure 
18a  are  the  observed  surface  temperatures  as  presented 
previously.  Immediately  apparent  is  an  upward  trend  in  the 
model  temperatures  relative  to  the  data,  resulting  in  large 
deviations  from  the  observations  by  the  9th  through  the  nth 
of  September.  In  fact,  least  squares  fits  of  the  data  and  the 
model  output  indicate  only  a  slight  upward  trend  in  the  data 
of  0 . 04 °C  over  the  11  day  period,  but  an  upward  trend  of 
2 . 93 °C  over  the  10.17  days  of  model  output.  This  corresponds 
to  an  average  increase  of  0.289°C  per  day.  This  trend  will  be 
discussed  in  the  following  chapter  with  respect  to  possible 
advective  effects. 

Also  of  note  in  Figure  18a  is  the  fact  that  while  the 
general  temporal  agreement  between  the  model  and  observed 
temperature  patterns  is  reasonably  good,  the  diurnal  cycle  is 
under-represented  by  the  Garwood  model.  Daily  increases  are 
sharp,  as  they  are  in  the  data,  but  in  all  cases  except  the 
9th,  exhibit  a  lower  magnitude  than  those  in  the  data.  Of 
greater  significance  in  the  overall  trend,  the  nocturnal 
cooling  in  the  model  is  much  less  than  observed,  especially  on 
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days  2,  3,  and  7  when  a  large  degree  of  actual  cooling 
occurred.  Also,  the  rate  at  which  the  cooling  occurred  in  the 
model  was  greatly  reduced  relative  to  that  seen  in  the 
observations.  Mixed  layer  cooling  on  the  9th  and  10th,  the 
days  with  heavy  cloud  cover,  was  more  in  line  with  that  seen 
in  the  data. 

Specific  characteristics  of  the  model  output  for  the 
purpose  of  comparison  with  surface  temperature  data  are 
offered  here.  Due  to  the  large  trend  in  the  data,  only  timing 
of  the  diurnal  events  and  magnitude  of  the  daily  temperature 
cycles  will  be  given.  The  Garwood  model  produced  mixed  layer 
temperature  minima  at  a  mean  time  of  1418Z  (0718  PDT)  and 
maxima  at  2206Z  (1506  PDT) .  Thus,  the  modeled  temperature 
minima  lag  the  data  by  an  average  of  44  minutes,  while  the 
peaks  actually  lead  those  in  the  data  by  an  average  of  57 
minutes.  The  day  to  night  temperature  falls  ranged  from 
0.14°C  to  l.01°C,  with  an  average  drop  of  0.42°C.  This 
compares  with  an  average  drop  of  0.96°C  in  the  data.  Thus, 
the  Garwood  model  cooled  the  mixed  layer  by  an  average  0.54°C 
less  than  that  seen  in  the  surface  temperature  data. 

Mixed  layer  depths  show  a  pattern  consistent  with 
these  thermal  results.  The  mixed  layer  can  be  seen  to  deepen 
the  most  on  the  days  when  the  cooling  is  the  strongest.  As 
previously  mentioned,  no  direct  measurements  of  mixed  layer 
depths  were  available.  However,  Figure  19  shows  the 
difference  between  the  observed  surface  temperatures  and  those 
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at  10  meters.  The  mixed  layer  can  be  seen  to  extend  to  at 
least  10  meters  depth  (at  least  briefly,  but  for  most  of  a  day 
on  more  than  one  occasion)  on  the  2nd  through  the  8th,  and 
again  on  the  10th.  By  contrast,  Figure  18b  indicates  that  the 
modeled  mixed  layer  deepens  below  10  meters  only  on  the  3rd, 
6  th,  7th,  and  8  th.  In  general,  then,  it  would  appear, 
consistent  with  the  lack  of  cooling  indicated  in  the  model 
temperatures,  that  this  model  also  somewhat  underestimates  the 
depth  to  which  the  diurnal  mixed  layer  deepening  occurs. 

Temperature  profiles  at  the  start  of  the  model  run  and 
at  specified  times  thereafter  are  presented  in  Figure  20. 
From  this  progression,  it  can  be  seen  that  the  stratification 
in  the  water  column  is  increasing  with  time,  such  that  by  192 
hours  (the  end  of  the  eighth  day  of  model  time)  significant 
warming  has  occurred  through  the  upper  13  meters.  These 
results  are  consistent  with  the  warming  trend  seen  in  the 
model  mixed  layer  temperature  results. 

2.  PWP 

.  The  PWP  model  results,  in  response  to  the  same  forcing 
as  used  in  the  Garwood  model,  with  friction  and  absorption  as 
defined  in  the  last  section,  are  shown  in  Figures  21a  and  b. 
A  surprisingly  similar  pattern  of  temperatures  and  mixed  layer 
depths  resulted.  Again,  there  is  an  upward  trend  in  the  model 
temperatures  relative  to  the  data,  although  somewhat  less  than 
that  seen  in  the  Garwood  model,  of  2.00°C  over  the  10.17  day 
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period  of  study,  or  0.196°C  per  day.  The  dampening  of  the 
diurnal  temperature  variation  is  even  more  extreme  in  this 
model,  with  very  little  cooling  occurring  on  any  nightly 
cycle,  and  with  too  little  diurnal  warming  as  well. 

PWP  produced  minimum  temperatures  at  a  mean  time  of 
1426Z  (0726  PDT) ,  quite  consistent  with  Garwood.  The  maximum 
also  agreed  very  closely  with  the  Garwood  results,  occurring 
at  an  average  2206Z  (1506  PDT) .  Thus,  there  is  again  a  lag  in 
the  model  temperature  troughs  of  an  average  52  minutes,  while 
the  daily  maxima  occur  an  average  60  minutes  earlier  in  this 
model  than  they  do  in  the  data.  The  even  more  inadequate 
nocturnal  cooling  in  this  model  is  evidenced  by  the  fact  that 
the  mean  day  to  night  mixed  layer  temperature  fall  was  only 
0 . 20°C.  The  range  of  temperature  drops  was  from  O.H8°C  on 
the  5th  to  0 . 313 °C  on  the  6th,  when  strong  winds  allowed 
maximum  shear  production.  Even  on  this  day,  however,  the 
model  fell  far  short  of  the  0.96°C  average  cooling  in  the 
data. 

As  would  be  expected  with  its  cooler  overall 
temperatures,  the  PWP  model  produced  mixed  layer  depths  that 
were  slightly  greater  than  those  produced  by  Garwood.  (Note 
that  the  PWP  model  outputs  mixed  layer  depth  in  integer  form 
in  accordance  with  the  vertical  bin  spacing  of  the  model, 
while  Garwood's  model  computes  h  precisely  from  the  TKE 
relationships.)  Following  the  strong  wind  day  of  the  6th,  in 
fact,  the  layer  deepened  to  16  meters  in  this  model.  The 
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mixed  layer  depth  reaches  10  meters  or  more  on  the  2nd,  4th, 
and  the  6th  through  the  8th.  This  would  seem  to  be  slightly 
more  in  agreement  with  the  observed  occasions  of  deepening  to 
10  meters  or  more. 

Again,  a  series  of  model  generated  temperature 
profiles  is  provided  in  Figure  22  to  give  a  more  complete  view 
of  the  model's  development  of  the  mixed  layer  with  time.  As 
in  the  Garwood  profiles,  stratification  can  be  seen  to  build 
with  time  such  that  by  the  end  of  the  eighth  day  warming  has 
occurred  down  to  16  meters.  In  PWP,  however,  the  profiles  are 
smoother  and  the  stratification  is  concentrated  between  15  and 
18  meters.  In  Garwood,  the  profiles  are  characterized  by 
significant  stratification  over  more  than  one  depth  band,  as 
can  be  seen  in  the  last  panel  of  Figure  20. 
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Figure  13a.  CTD  Temperature  Profiles  From  9  September  at 
Various  Points  Within  Monterey  Bay.  Profile  #2  was  used 
for  model  runs,  as  discussed  in  text. 
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Figure  13b.  Initial  Temperature  Profile  Used  by  the 
Models . 
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Figure  14.  PWP  Mixed  Layer  Temperatures  With  Varying 
Friction. 
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September,  1992 


Figure  15.  Step  wind  Function  (Top)  Applied  to  Garwood 
Unsteady  Model  and  Resulting  Total  TKE  (Bottom)  Versus 
Time. 
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Figure  16b.  PWP  Mixed  Layer  Temperatures  With  and  Without 
Clouds . 
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September,  1992 


Figure  17.  Total  Absorption  Profiles  Computed  by  Both 
Models.  Garwood  Absorption  is  Dashed  and  PWP  is  Solid. 
Note  that  PWP  Breaks  the  Absorption  into  Red  and  Blue 
Profiles  in  the  Ratio  0.6  Red  and  0.4  Blue. 
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Figure  18a.  Final  Garwood  Mixed  Layer  Temperatures  Versus 
Observed  Sea  Surface  Temperatures. 
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Figure  18b.  Garwood  Mixed  Layer  Depths. 
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September,  1992 


Figure  19 .  Observed  Sea  Surface  Temperatures  Minus  10 
Meter  Temperatures.  Values  Near  Zero  Indicate  Periods  of 


Deepening  to  Greater  than  10  Meters. 
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September,  1992 


Indicated  at  the  Top  of  Each  Profile. 
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Figure  21a.  Final  PWP  Mixed  Layer  Temperatures  Versus 
Observed  Sea  Surface  Temperatures. 
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September,  1992 


Figure  21b.  PWP  Mixed  Layer  Depths.  Note  That  Depths  are 
in  Whole  Numbers,  Since  They  are  Inferred  from  a  Density 
Profile  Which  is  Subject  to  the  Model's  Vertical  Bin 
Spacing. 
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September,  1992 


Figure  22.  Sequence  of 
Model . 
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Profiles  from  the  PWP 


VI.  DISCUSSION 


A.  RELATING  MODEL  TRENDS  TO  ADVECTION 

As  stated  in  the  last  chapter,  both  models  exhibited  an 
upward  trend  in  mixed  layer  temperatures  relative  to  the 
observations.  These  trends  were  0.289  and  0.196°C  per  day  for 
Garwood  and  PWP,  as  presented  previously.  The  least  squares 
fits  used  to  compute  these  trends  are  shown  in  Figures  23a  and 
b.  If  it  is  assumed  that  this  trend  is  due  entirely  to  the 
fact  that  no  advective  effects  are  included  in  the  models  and 
that  any  advection  occurring  is  in  the  vertical,  then  an 
estimate  of  the  magnitude  of  vertical  advection  can  be  made, 
according  to  the  relationship 


dT 

at 


(37) 


Using  a  vertical  temperature  gradient  from  the  initial  profile 
of  -2.4°C  over  the  top  21  meters  (the  depth  region  with 
significant  stratification),  or  -0.114°C/m,  produces  vertical 
velocities  of  2.53  and  1.72  meters  per  day  for  Garwood  and 
PWP,  respectively.  These  values  are  within  the  range  of 
upward  velocities  expected  in  an  upwelling  region  (Huyer, 
1983)  .  Therefore,  it  is  reasonable  to  suggest  that  the  upward 
trend  in  both  models'  results  could  be  entirely  due  to  the 
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fact  that  no  advection  was  included  in  the  versions  of  the 
models  used  here. 

Figure  24  shows  the  least  squares  fit  lines  of  the  raw 
data  for  surface,  10,  20,  and  40  meter  temperatures.  A  clear 
downward  trend  in  the  subsurface  data  suggests  that  colder 
water  was  being  advected  into  the  upper  40  meters,  which  could 
have  contributed  to  maintaining  fairly  constant  surface 
temperatures  despite  forcing  that  should  produce  a  warming 
trend.  This  is  another  indication  that  vertical  and/or 
horizontal  advection  was  significant  during  the  study  period. 

B .  DIURNAL  CYCLE 

Apart  from  the  trends  seen  in  the  models  relative  to  the 
surface  temperature  data,  the  ^ther  major  deviation  of  the 
models,  as  mentioned  in  the  last  chapter,  is  in  the  magnitude 
and  characteristic  shape  of  the  diurnal  cycle.  This  can  be 
seen  clearly  if  the  aforementioned  trends  are  removed  from  the 
model  results,  as  done  in  Figures  25a  and  b.  The  models  do 
produce  a  sharp  increase  in  temperature  each  morning  as  the 
sun  comes  up,  quite  similar  to  that  seen  in  the  data,  although 
the  magnitude  of  this  rise  is  in  most  cases  much  less  than 
that  seen  in  the  data.  After  the  models  reach  their  peak  sea 
surface  temperature,  which  occurs  at  approximately  1506  PDT 
(one  hour  prior  to  the  average  time  in  the  data) ,  they  begin 
a  decline  in  temperature  which  is  much  more  gradual  than  the 
observations.  This  is  true  in  both  models,  although  more  so 
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in  PWP .  This  reduced  downward  slope  in  mixed  layer 
temperatures  in  the  Garwood  case  results  in  nightly  decreases 
of  an  average  0.54°C  less  than  in  the  data.  The  shortfall  in 
nocturnal  cooling  in  PWP  is  even  greater,  at  0.76°C. 

This  discrepancy  between  the  models  and  observations 
dictated  further  analysis  of  model  behavior  in  response  to 
varying  forcing.  In  an  effort  to  determine  how  the  models 
were  responding  to  heating  and  wind  stress  inputs 
independently,  two  additional  test  cases  were  run.  The 
results  of  the  model  runs  presented  in  the  previous  chapter 
and  detrended  in  Figures  25a  and  b  show  how  the  models 
reproduced  mixed  layer  temperature  in  response  to  the  combined 
effects  of  varying  wind  stress  and  heat  flux.  In  the  first 
test  case  the  wind  stress  is  held  constant,  while  the  heat 
flux  undergoes  the  usual  diurnal  variation.  This  is  similar 
to  what  has  been  done  in  open  ocean  studies.  Secondly,  in  an 
effort  to  measure  the  model  response  to  diurnal  variations  of 
wind  stress  alone,  winds  were  allowed  to  vary  according  to  the 
observations  while  the  total  heat  flux  was  held  constant. 

1.  Constant  Wind  Forcing 

The  constant  wind  cases  were  run  using  constant 
eastward  and  northward  winds  of  3.672  and  -1.038  m/s, 
respectively.  These  were  the  mean  values  observed  over  the  11 
day  period,  corresponding  to  a  mean  wind  toward  117°T,  as 
mentioned  in  Chapter  IV.  Temperature  results  from  both  models 
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are  presented  in  Figure  26a,  while  mixed. layer  depths  are 
shown  in  Figure  26b.  An  extremely  regular  pattern  of  mixed 
layer  temperatures  is  apparent,  as  would  be  expected  from  the 
relatively  regular  solar  heating  pattern.  Mixed  layer  depths 
remain  shallow  in  the  absence  of  afternoon  wind  peaks.  The 
Garwood  model  exhibits  significantly  more  diurnal  variability 
than  PWP,  warming  an  average  l.44°C  each  day  as  compared  to 
only  0.67°C  per  day  in  PWP.  Both  models  cool  the  mixed  layer 
each  night  by  an  amount  less  than  they  warmed  it  during  the 
day.  In  the  case  of  the  Garwood  model,  this  cooling  averages 
0.91°C  per  day,  implying  a  net  gain  of  mixed  layer  temperature 
each  day  of  0.53°C.  PWP  cools  by  an  average  0.18°C  per  day, 
resulting  in  a  similar  gain  of  0.49°C.  These  substantial 
daily  heat  gains  when  the  model  is  forced  with  constant  winds 
of  this  magnitude  produce  the  large  upward  trends  seen  in 
Figure  26a.  This  suggests  that  daily  intensifying  winds  are 
an  essential  element  in  predicting  the  thermal  structure  in  a 
sea  breeze  influenced  region.  Without  the  increased 
generation  of  turbulence  from  wind  stress  that  occurs  in  the 
afternoons  and  early  evenings,  the  large  downward  heat  flux 
which  is  occurring  at  those  times  creates  an  unrealistically 
shallow  mixed  layer,  and  thus  excessive  mixed  layer  warming. 
Clearly,  the  timing  of  the  peak  in  wind  stress  is  critical. 
The  fact  that  it  occurs  at  a  time  when  solar  insolation  is 
high  allows  the  incoming  heat  to  be  distributed  over  a  deeper 
layer  than  would  otherwise  be  the  case.  If  winds  were  light 
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in  the  afternoons  and  intensified  at  dusk,  for  example,  much 
greater  daily  warming  and  nightly  cooling  would  occur. 

Of  note  in  the  constant  wind  case  is  the  fact  that 
temperature  peaks  occur  at  an  average  time  of  0014Z  and  0057Z 
(1714  and  1757  PDT)  in  Garwood  and  PWP,  respectively.  Recall 
that  the  observed  surface  temperatures  peaked  at  an  average 
time  of  1603  PDT,  while  the  model  runs  forced  with  varying 
wind  stress  and  heat  flux  peaked  at  an  average  time  of  1506 
PDT.  Thus,  in  the  constant  wind  case,  there  is  a  lag  with 
respect  to  the  data,  rather  than  the  lead  that  was  seen  in  the 
total  forcing  cases.  Allowing  the  wind  to  vary,  then, 
resulted  in  a  phase  shift  of  2  to  3  hours,  depending  on  the 
model,  since  increased  afternoon  mixing  was  able  to  begin 
cooling  the  layer  sooner  than  when  the  winds  were  held 
constant. 

2 .  Constant  Heat  Flux 

For  these  cases,  the  wind  stress  computed  from 
observed  eastward  and  northward  wind  components  as  discussed 
previously  was  used,  along  with  a  constant  downward  surface 
heat  flux  of  157.7  W/m2,  to  force  the  models.  Again,  this  was 
the  mean  value  over  the  study  period  and  is  composed  of  a  mean 
solar  insolation  of  231.3  W/m2  (including  nighttime  periods) 
and  a  mean  heat  loss  of  73.6  W/m2.  Mixed  layer  temperature 
and  depth  results  for  these  runs  for  the  two  models  are 
presented  in  Figures  27a  and  b.  These  results  for  the  Garwood 
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model  very  nearly  duplicate  those  that  resulted  from  the 
complete  forcing  case  presented  in  the  last  chapter,  with 
slight  deviations  occurring  only  over  the  last  half  of  the 


period.  PWP  also  produced  a  thermal  pattern  which  more 
closely  resembled  that  of  the  last  chapter  than  did  the 
constant  wind  stress  case.  However,  in  PWP,  the  diurnal 
variability  is  even  more  damped  out  than  it  was  previously. 

Particularly  with  respect  to  the  Garwood  model,  these 
observations  suggest  that  the  details  of  the  wind  stress  used 
to  force  the  models  is  a  greater  determinant  of  the  results 
than  the  heat  flux.  The  fact  that  the  total  forcing  case  and 
the  case  with  constant  downward  heat  flux  produced  nearly 
identical  mixed  layer  temperature  patterns  is  quite 
surprising,  and  might  suggest  that  the  model  does  not  respond 
to  the  diurnal  heat  flux  cycle  in  a  manner  similar  to  the  real 
ocean. 

C.  POSSIBLE  REASONS  FOR  MODELS'  UNDERESTIMATES  OF  NOCTURNAL 

COOLING 

Mechanisms  that  might  cause  the  models  to  underpredict  the 
large,  sharp  temperature  falls  that  are  seen  on  most  nights  in 
the  data  are  discussed  next. 

X.  Penetrative  Convection 

After  the  sun  sets,  surface  heat  flux  from  the  ocean 
to  the  atmosphere  generates  higher  density  water  and  thus 
convective  instability.  Both  models  deal  with  this  by  mixing 
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that  water  downward  until  static  stability  is  returned 
throughout  the  water  column.  One  suggestion  to  explain  the 
reduced  cooling  in  the  PWP  model  relative  to  the  observations 
each  night  is  that  the  cool  water  at  the  surface  actually 
descends  in  plumes  which  have  associated  vertical  momentum. 
This  causes  them  to  overshoot  the  base  of  the  mixed  layer, 
entraining  a  larger  quantity  of  dense  water  up  into  the  layer. 
This  would  result  in  a  cooler  mixed  layer  by  morning  than 
would  occur  without  this  process,  and  would  presumably  cause 
the  mixed  layer  to  deepen  more  rapidly  at  the  onset  of  upward 
buoyancy  flux.  PWP  does  not  include  this  penetrative 
convection  mechanism.  Garwood  does  include  such  a  mechanism 
and  this  could  explain  its  increased  nocturnal  cooling 
relative  to  PWP.  This  mechanism  is  suggested  by  Large,  et  al. 
(1993)  . 

2.  Steady  State  TKE  Assumption 

Garwood's  model  assumes  that  the  TKE  is  balanced  at 
every  time  step.  As  discussed  previously,  an  unsteady  model 
(which  handles  only  deepening  of  the  mixed  layer)  was  run  with 
step  function  wind  forcing.  As  seen  in  Figure  15,  the  TKE 
reaches  a  steady  state  in  about  an  hour,  but  during  that  hour 
TKE  exceeds  the  value  to  which  it  assymptotes.  The  presence 
of  these  short  term  transients  in  the  ocean,  which  could  be 
generated  by  constantly  varying  buoyancy  forcing,  as  well  as 
wind  stress,  could  contribute  to  mixing  that  is  not  reproduced 
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by  the  steady- state  model.  This  additional  mixing,  then, 
would  result  in  cooler  mixed  layer  temperatures. 

3.  Model  Stratification 

At  initialization,  it  is  known  chat  a  fairly  accurate 
temperature  profile  is  being  used  by  the  model.  However,  as 
the  model  mixed  layer  temperature  increases,  while  the  profile 
below  the  layer  remains  the  same,  there  is  a  corresponding 
increase  in  the  degree  of  stratification  in  the  model  as  time 
progresses.  This  excess  density  contrast  makes  subsequent 
mixing  more  difficult.  That  is,  a  greater  amount  of 
turbulence  is  necessary  for  the  same  amount  of  entrainment. 
A  visual  representation  of  this  temperature  profile  change  was 
seen  in  Figures  20  and  22.  The  fact  that  on  the  first  night 
of  the  period,  before  this  effect  could  build,  the  models 
produced  a  temperature  fall  in  line  with  the  data  lends 
support  to  this  idea  as  a  potential  contributing  factor. 

4 .  Diurnal  Advection 

Another  postulated  mechanism  for  producing  the  large 
and  sharp  temperature  decreases  seen  in  the  data  at  night  is 
a  diurnally  varying  advection.  Although  currents  were  not 
specifically  studied  here,  Foster  (1993)  demonstrates  from  HF 
radar  measurements  of  surface  currents  in  Monterey  Bay  that 
the  spectrum  of  surface  currents  is  dominated  by  the  diurnal 
period,  presumably  driven  by  the  sea  breeze.  It  is  possible 
that,  in  response  to  the  very  regular  variation  of  the  surface 
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currents,  there  could  be  a  similarly  regular  pattern  of 
horizontal  (or  vertical)  advection  occurring  within  the  bay. 
In  order  for  horizontal  advection  to  produce  large  temperature 
changes  at  Ml,  there  must  be  a  significant  horizontal 
temperature  gradient  in  the  area,  which  is  being  forced  across 
the  mooring  site  on  a  diurnal  time  scale.  Given  the  presence 
of  upwelled  waters  in  the  region  and  the  location  of  Ml  near 
the  mouth  of  the  bay,  it  seems  possible  that  such  diurnally 
varying  advective  effects  could  be  occurring.  Without 
detailed  surface  temperature  data  at  high  temporal  resolution 
in  the  area,  it  is  not  possible  to  test  this  hypothesis. 
Diurnal  variations  in  vertical  advection  would  require 
horizontal  divergence/convergence  of  the  diurnal  surface 
currents.  This  hypothesis  also  requires  additional  data  to 
test . 
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Figure  23a.  Least  Squares  Fit  to  Garwood  Temperatures 


September,  1992 


,  1992 


Surface 
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September,  1992 


Figure  25a.  Detrended  Garwood  Mixed  Layer  Temperatures 
(Solid)  Versus  Observed  Sea  Surface  Temperature.  Note 
Different  Temperature  Scale  from  Previous  Plots. 
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September,  1992 


Figure  25b.  Detrended  PWP  Mixed  Layer  Temperatures  Versus 
Observed  Sea  Surface  Temperatures. 
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,  1992 


Figure  26a.  Mixed  Layer  Temperatures  with  Constant  Wind, 
Varying  Heat  Flux. 
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September,  1992 


Figure  27a.  Mixed  Layer  Temperatures  with  Constant  Heat 
Flux,  Varying  Wind.  Note  Similarity  of  Garwood  Pattern  to 
Previous  Result. 
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September,  1992 


Figure  27b.  Mixed  Layer  Depths  with  Constant  Heat  Flux, 
Varying  Wind. 
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September,  1992 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  MIXED  LAYER  BEHAVIOR  IN  MONTEREY  BAY 

The  behavior  of  Che  oceanic  mixed  layer  in  a  region 
influenced  by  diurnally  intensifying  winds  is  studied.  The 
study  area  is  Monterey  Bay,  California,  in  which  a  sea  breeze 
is  clearly  shown  to  be  strongly  influencing  the  temperature 
and  depth  of  the  mixed  layer  during  September,  1992.  This 
diurnal  wind  stress  variability  creates  an  interesting 
interaction  with  the  diurnal  heat  flux  present  in  other  mixed 
layer  studies.  The  wind  stress  reaches  its  peak  at  an  average 
time  of  1730  PDT.  The  downward  heat  flux  peaks  at  about  1300 
PDT.  The  sea  surface  temperature  data  presented  in  Chapter  IV 
shows  a  pattern  of  large  day  to  night  temperature  swings,  with 
the  peaks  occurring  at  approximately  1600  PDT.  In  the 
evening,  decreasing  downward  heat  flux  and  large  wind  stress 
produce  cooling,  which  is  underpredicted  by  the  models. 

The  models'  mixed  layer  de  th  is  also  shown  to  undergo  a 
diurnal  variation  influenced  by  the  sea  breeze.  It  is  seen  to 
deepen  quickly  in  the  evenings  when  winds  are  still  quite 
strong  and  reach  its  greatest  depth  just  prior  to  sunrise 
after  a  night  of  convective  overturning.  Shallowing  in  the 
morning  occurs  very  suddenly  in  the  presence  of  the  light 
winds  seen  at  those  hours. 
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B.  ONE -DIMENSIONAL  MIXED  LAYER  MODELS  IN  SEA  BREEZE  REGIONS 

The  one -dimensional  mixed  layer  models  of  Garwood  (1977) 
and  Price  et  al.  (1986)  are  used  to  evaluate  the  capabilities 
of  such  models  in  this  coastal  environment.  The  models  do  a 
good  ]ob  of  predicting  the  phase  of  the  mixed  layer 
temperature  pattern  seen  in  the  data.  Compared  to 
observations,  a  lead  of  about  one  hour  is  seen  in  the  model 
temperature  peaks,  with  a  lag  of  less  than  an  hour  in  the 
minima.  The  daily  increases  in  mixed  layer  temperature  are 
very  similar  in  slope,  although  frequently  with  reduced 
magnitude,  in  the  models  to  what  is  observed  in  the  data. 

Both  models  exhibit  two  major  differences  from  what  is 
seen  in  the  observations  over  this  period.  The  first  of  these 
is  the  upward  trend  in  mixed  layer  temperatures  relative  to 
the  data  (0.29°C  per  day  in  Garwood  and  0.20°C  per  day  in 
PWP) .  This  is  postulated  to  be  primarily  the  result  of  the 
fact  that  no  advection  is  included  in  either  model . 
Reasonable  values  of  vertical  advection  are  obtained  using  the 
assumption  that  all  of  the  missing  advection  is  in  the 
vertical.  The  other  characteristic  of  the  model  results  that 
differs  from  the  data  is  the  fact  that  the  cooling  at  night  is 
reduced  in  both  magnitude  and  rate.  This  is  true  in  both 
models,  although  the  Garwood  model  produces  a  better  diurnal 
cycle  than  does  PWP.  This  lack  of  cooling  is  seen  to  be  true 
to  the  same  degree  when  the  Garwood  model  is  forced  with  a 
constant  downward  heat  flux.  These  results  suggest  that,  at 
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least  in  the  presence  of  diurnal  wind  stress  variation,  the 
upward  heat  flux  phase  of  the  daily  buoyancy  forcing  cycle  has 
little  influence  on  the  model.  Two  mechanisms  which  could 
produce  the  large,  sharp  temperature  drops  observed  each  night 
in  the  data  are  penetrative  convection  and  diurnally  varying 
advection.  The  former  is  not  present  in  PWP  and  the  latter  is 
not  included  in  either  model.  Transients  in  the  turbulent 
kinetic  energy  produced  by  the  constantly  varying  forcing  are 
not  included  in  the  Garwood  model,  which  assumes  steady  state 
TKE  at  each  time  step.  This  and  the  increasing  stratification 
in  the  models  as  time  progresses  would  reduce  the  amount  of 
mixing  and  cooling  that  occur  in  the  models  relative  to  that 
in  the  ocean. 

Overall,  it  is  concluded  that  the  turbulence  budget  type 
of  model  reflects  the  real  ocean  in  this  environment  to  a 
greater  degree  than  does  the  Richardson  number  instability 
type,  based  on  the  models  studied  here.  However,  both  types 
apparently  suffer  from  their  one-dimensionality  in  a  coastal 
region  where  vertical  and  horizontal  advective  effects  are 
apparently  present. 

C.  RECOMMENDATIONS 

Since  it  is  possible  to  add  vertical  advection  to  a  one¬ 
dimensional  model,  it  is  recommended  that  this  be  done  in  an 
area  such  as  this  where  these  effects  are  clearly  not 
negligible.  By  adding  the  advection  at  each  model  time  step, 
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the  stratification  would  not  build  as  much  as  was  seen  in  this 
study  ana  a  more  realistic  trend  would  result.  Also,  it  is 
possible  that  this  would  also  allow  a  more  realistic 
prediction  of  nocturnal  cooling,  particularly  later  in  the 
model  run. 

In  order  to  investigate  the  possibility  of  additional 
mixing  due  to  transients  of  TKE ,  it  would  be  beneficial  to  run 
a  non- steady  state  model  at  a  shorter  time  step  (say,  10 
minutes)  for  the  entire  period.  It  would  be  necessary  to  use 
a  model  which  allows  both  deepening  and  shoaling  of  the  mixed 
layer. 

This  study  brings  up  a  variety  of  possible  future  efforts 
related  to  coastal  mixed  layer  physics.  Conducting  a  similar 
study  with  this  data  during  another  time  of  year  (January,  for 
example)  when  the  sea  breeze  is  not  well  established  would 
provide  an  interesting  comparison.  Since  upwelling  is  reduced 
or  absent  along  the  coast  during  this  time,  and  the 
temperature  field  is  more  uniform  in  the  horizontal,  it  is 
possible  that  advective  effects  would  also  be  reduced,  and  the 
models  would  more  accurately  reflect  the  data. 

Use  of  the  surface  current  data  from  high  frequency  radar, 
such  as  that  used  by  Foster  (1993)  ,  would  allow  convergence  at 
Ml  to  be  computed,  thus  making  possible  a  determination  of  the 
vertical  velocity  near  the  surface.  This,  together  with  the 
vertical  temperature  gradient,  can  be  used  to  calculate 
vertical  advection.  To  determine  the  variability  of 
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horizontal  advection,  it  would  be  useful  tc  obtain  an  improved 
picture  of  sea  surface  temperatures  within  the  bay.  This 
could  be  done  through  a  combination  of  extensive  in  situ 
measurements  and  the  use  of  satellite  derived  sea  surface 
temperatures.  Combining  this  information  with  velocity 
results  produced  by  the  models  may  provide  further  insights 
into  the  three-dimensional  forcing  occurring  at  Ml  and  other 
coastal  locations. 

Finally,  a  more  complete  picture  of  mixed  layer  behavior 
should  be  obtained  through  a  series  of  upper  ocean  CTD  or 
bathythermograph  casts  collocated  with  continuous 
meteorological  observations.  These  casts  should  be  of 
sufficient  frequency  to  allow  resolution  of  changes  on  an 
hourly  or  smaller  time  scale.  This  would  eliminate  the 
problem  encountered  in  this  study  in  which  temperature  data 
was  available  from  only  the  surface,  10  and  20  meters.  This 
would  also  allow  a  more  accurate  description  of  density 
variations,  since  salinity  would  be  included. 
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